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ABSTRACT 


Seventy-five measurements of the differential elastic 
electron-deuteron scattering cross section were made for 
values of ae ranging trem 0.05 to 0.35 fm ° at the NPS LINAC 
with electron scattering angles of 60° to 120°. The charge 


radius of the deuteron was determined to be te = 1.9641 + 
d 
0.0074 fm. Values of the neutron charge form factors were 
calculated using Lomon-Feshbach Models 1, 5, and 15 and 
relativistic corrections to the deuteron wave functions. 
Lomon-Feshbach Model 15 gave —*5 G, (0) = 0.0200 + 0.0058 fm’, 
dq n 

in closest agreement with the value of Go (0) determined 

dg n 
by thermal neutron work. A Rosenbluth plot of the data at 


5c = 0.2 shows good agreement with theory. 
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I. INTRODUCTION 


Pee OUTLINE OF THE EXPEREMENET AND ENTRODUCTION OF FORM FACTORS 
The primary purpose of this work was to measure the charge 
form factor of the neutron, GE (q7), at low values of q? (the 
square of the Lorentz eee four-momentum transfer from 
fe scattered electron to the neutron) in order to resolve 
Rai ehces among several earlier experiments [1,2,3,4]. 


Secondly, by comparing the Gr obtained with the use of 


several different deuteron wave Foner ene to the Gp jace= 
dicted by accurate thermal neutron experiments, eee select 
the S state and D state combinations of the deuteron wave 
functions that are in best agreement with experiment. ~ 

Thirdly, measurements of the deuteron cross section at 
low values of q? permit the rms charge radius of the deuteron 
to be determined accurately. 

These three objectives were achieved by measuring the 
elastic electron - deuteron differential scattering cross 
section. The elastic electron - proton cross section was 
measured in order to normalize the measured deuteron cross 
section. The electron beam was produced by the Naval 
Bi dradivate School linear accelerator (LENAC). 

Since the primary purpose was the measurement of the charge 
form factor, it is useful to consider briefly here what a form 
factor is. In Section II a more detailed explanation will be 
given. The cross Seetiion, 22, for a relativistic electron 


an 


scattering from a nucleus that has zero spin and magnetic 


10 





moment and consists only of a point charge has been calcu- 
lated by Mott [5]. The experimental cross section is re- 


lated to the so-called Mott cross section by 


ee ee ~ = jae 6" (q") (I-1) 
Wmere G is the form factor. If the nucleus has no structure, 
iie., tt is a point, G = 1 at all values of aoe Ef it has 
memucture, G = 1 for q? equal to zero and G is less than one 
fOr a? Greater than zero. <A reason why Gta) is called the 
nom fL£actor [3] is seen when the cross section for elastic 
Bewecering of an electron from a nucleus of charge Ze with 


finite extent is written [6] as 


(eee Com [(eurerS Fare He ae 
Gin: oA ee een. 7 MO tt Ze d 


where P(r) is the nuclear charge distribution and chq is 
the three-momentum transferred to the nucleus. The term in 
brackets is the form factor G and is the Fourier transform 
of the charge distribution. If one calculates the theo- 
retical cross section for elastic scattering, it is seen 
Ew@at the charge distribution involves the preduct of ae 
initial wave function of the scatterer at rest in the labo- 
ratory frame ne the final wave function of the scatterer 
after it has absorbed the recoil momentum chq. In Section 
II it will be shown that a relativistic correction to the 
deuteron wave function is required because of its recoil 
velocity. This is a small (a few percent) but important 


Gerrection in this work. 


del 





L£ we consider elastic scattesing of an eclleeteons iron 
an infinitely massive (that is, a nucleus which does not 
recoil) nucleus of charge Ze Equation I-2 will correctly 
describe the scattering; P(r) is thensthe static charge 
amet ribution. usecomeeee light nuclei where recoil is im- 
portant, such as deuterium and hydrogen which were used as 
targets in this work, the intuitive concept of the form 
factor as the Fourier transform of the charge distribution 
is not unambiguous. For example, if the charge distribution 
Ort the proton were spherical in one frame of reference then 
in another frame the Lorentz contraction of the distribution 
would cause it to be nonspherical. To describe form factors 
in a way more useful in theoretical work, they are considered 
as invariant functions of the four-momentum gq rather than the 
three-momen tum cha, where are = heeee ~ (E,-E,)* and Es 
and E, are the initial and final energies. When the recoil 
is small then q? eS a7 andthe conmeecpt of the form factor as 
the Fourier transform of the charge distribution has useful- 
ieee Asewill beeseen in Sectzon II the inwaryant form 
factors have other theoretical meanings in terms of strong 
interaction dynamics. 

Mhere is another complication. Particles and nuclel may 
have magnetic moment distributions as well as charge distri- 
butions. These distributions may or may not have the same 
form as the charge distributions. Thus, a particle or nuc- 
leus that has a magnetic moment (but no quadrupole or higher 
moments) will have to be described by two independent form 


factors. The decomposition into electric and magnetic 


12 





distributions depends upon the frame of reference. For ex- 
ample, the elastic scattering of an electron from a massive 
nucleus with charge Ze and no magnetic moment will be cor- 
rectly described in the laboratory frame by Equation I-2 if 
the nucleus is at rest in the laboratory frame. If, however, 
the calculation is performed in the electron's frame of refer- 
ence then the moving nucleus will also have a magnetic field 
and thus an apparent magnetic moment. In this work we are 
ultimately interested in the electric form factor of the 


neutron. A particularly useful frame, the Breit or "brickwall 


frame", separates the electromagnetic interaction of the nuc- 
leons into two form factors, Gas the electric form Eactor, 
and G,, the magnetic form factor. Sachs et al [7,8] showed 


that these are particularly useful form factors since Gh and 


Gy can be expressed as the Fourier transforms of the spatial 


distributions of charge and magnetization in the Breit frame. 


Wieeereit frame is the frame in which ae = aa 7 lve eeieae 


has ES = Ee and ‘that’ has the “absolute valiiG@ of,the momentum 


of the incoming nucleon the same as the absolute value of the 
momentum of the outgoing nucleon and oppositely directed [9]. 
Consider the electric form factor of a nucleon in the 


Breit frame, en 
‘ ig°r 43) 
e a 


etieg) = \ p(x) 2 


which for a spherically symmetric particle may be written as 


eo 


\ Sin (45) Pr spiG@aiar: 
O 


aot 


2 
GC, (q™) aa 
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For small values of ne the sine may be expanded to give: 


one 


1 


where ee is the root mean square radius of the charge 
distribution. A similar expansion may be made for the mag- 
netic distribution to obtain the magnetic oe Thesproken 
and neutron each have different form factors. Thus, four 


form factors are required to describe them: Gr ; Gy ’ Guy ' 
n 


and G, . The eat obtained by using Equation I-3 is deter- 

n 
mined in the Breit frame, but since Equation I-3 is only good 
at low values of mee the ane value determined is generally 


quoted as the nucleon'’s charge radius. 


Pee etnOD OF OBTAINING THE NEUTRON ELECTRIC FORM FACTOR, Gp 
N 


The process of obtaining G from the measured elastic 


L 
n 


electron-deuteron cross section will be outlined here and 
discussed in detail in Section II. First, the elastic elec- 
tron-proton cross section was measured at a given scattering 
angle and scattered energy using hydrogen gas as a target. 
The hydrogen was then replaced with deuterium gas and the 
elastic electron-deuteron cross section was measured at the 
same scattering angle and with the same scattered energy. 
The incident electron energy was adjusted to give this con- 
dition. All geometries, beam currents, gas pressures, and 
operating procedures were kept as constant as possible while 
both the proton and the deuteron data were taken. In fact, 
the LINAC was operated continuously from the time an experi- 


ment started until it was completed with both cross sections 


14 : 





measured. The time negQuired for Such “a“"runm 9 vari1ca. trom 
about 15 to 40 hours. 

Because absolute cross sections cannot be measured here 
owing to systematic uncertainties, the ratio of the cross 
sections, ot oP, was formed. In this ratio, many systematic 
uncertainties cancel since the cross sections are measured 
under the same conditions. The absolute proton cross section 
has been measured in previous experiments [10,11,12], and 
these data have been fit by deVries [10] to formulas from 
which the cross section can be calculated. When the experi- 
mental ratio above is multiplied by the absolute proton cross 
section, the absolute deuteron cross section is obtained. 


Several steps are required to extract Gh £rom this cross 
n 


section. The experimental cross section is first divided by 
the Mott cross section of the deuteron to obtain the deuteron 
fou Lactor, Go Gg accounts for all the structure an the 
deuteron, i.e., the electric and magnetic distributions 
within the neutron and proton as well as the convection of 
charge due to their motion within the deuteron and the mag- 
netic moment caused by this convection of charge. Following 
Jankus [13], the deuteron was treated as though its two con- 
stituent particles, the neutron and proton, were point par- 
ticles. The distribution of charge predicted by the deuteron 
wave function is then contained in a term called the deuteron 
charge structure factor, which depends only on a and the wave 
function of the deuteron ground state. E.L. Lomon [14] has 
provided three selected sets of values, from his larger group 


of sets of values, for these deuteron charge structure 
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factors. Each of these is based on a different deuteron wave 
function obtained from ntieleon=nweleon dawa [15], 

Drickey and Hand [3] experimentally verified a proportion- 
elity between the magnetic €orm factor ef the @ieWteroemmand its 
Garge form factor, called a “scaling law". That relationship 
was used to remove that part of the deuteron form factor 
caused by the intrinsic magnetic moments of the proton and the 
neutron as well as that part of the deuteron magnetic moment 
caused by the convection of charge. 

mt this point only the extended charge distribution of 
the proton has not been accounted for in the extraction pro- 
cess. This distribution is accounted for by the proton elec- 
tric form factor derived from deVries' data (see Appendix A). 


Mae wmeutron electric form factor, Gr , has thus been extracted 
n 


from the deuteron cross section. 

The above procedure is justified since the deuteron is a 
weakly bound structure, i.e., the proton and the neutron spend 
most of the time outside the range of nuclear forces and one 
can adequately describe the charge and current distribution 
Wathin each nucleon as if they were free [13]. 

Early in this work, before the gas target system was 


operating, efforts were made to measure Gp by a process simi- 
n 


lar to that described above but using polyethylene targets of 
CH, and cD. to obtain the cross sections. There were several 
problems with that work but the most important limitation was 
the scattering from carbon which complicated the data reduc-~ 


tion. The solid target data which were obtained were not as 


reproducible as the gas target data and are not included in 


this work. 
16 





€. RESULTS OF PREVIOUS EXPERIMENTS AND OCF wilco vere. 


in this. work thesneutronvchangestocrmetacto uaa was 


E é 
n 
measured for values of ne ranging from 0.05 .to O23 oe Z 


taeeslope of o versus q? was obtained by a linear least 


squares fit to ee data. 

The slope in the limit that q?o O has been measured in 
many experiments [16,17,18] which measured the interaction of 
thermal neutrons with atomic electrons. The most accurate 


emo most recent measurement of this kind is that of Krohn 


and Ringo [4], who obtained 


2, Gc, (0) = 0.0193 + 0.0004 £m*. (1-4) 
nN 


dq E 
This is probably that best determined of all nucleon form 


factor parameters [19]. Drickey and Hand [3] measured Gp 
n 


by measuring the ratio of the elastic electron~deuteron 
cross section to the elastic electron~proton cross section 


using liquid hydrogen and deuterium for targets and found 


that Gp = 0200 +,0.01 over a rangesot q? fremme0. Syren Ze2 ie 
n 


That is, they found 


G6, ~0, 

dg n 
Wateh iS in see eae with the thermal neutron measurements. 

Casper and Gross [20] developed relativistic corrections 
to electron-deuteron scattering and reanalyzed the data of 
Drickey and Hand. Using deuteron wave functions derived from 
the Feshbach~Lomon nucleon-nucleon model they found that the | 
data of Drickey and Hand indicated 
d 2 


Bez Cn (0) = 0.0153 + 0.0067 fm*, 
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which was consistent with the thermal neutron slope. 
Bumiller, et al. [21,2] extended the elastic electron-deuteron 
scattering data to a o£, 0.107 = q? = 0.8 fm, using CH., CD., 
and C targets and substantiated the conclusions of Casper and 
Gross byeobtaining) a slopesct,0.0179 +30 003G Unc igiekeas 
emmor of + 0.0036 is incorrect. It should have been quoted as 
+ 0.0089. Using the same equipment as Bumiller, Mader [1] 
continued the elastic electron-deuteron scattering work but 
used a gas target to avoid the carbon subtraction procedures 
which were necessary with the solid targets used by Bumiller 
eCtval. Although Mader‘'s data are as accurate as Refs. 21 and 
2, his measurements showed that i Gp ~*~ 0 over a range of 

n 
0.05 § q? = 0.60 an in this ae we have resolved this 


discrepancy by measuring the G slope to be 


E 
n 


2. G, (0) = 0.0200 + 0.0058 £m”, 


dq n 
witch then compares favorably to the value of 0.0193 + 0.0004 
fim" as measured directly by the neutron-electron scattering 


experiments. 


is 





Er. Vlaiso 


AS ELASTIC ELECTRON-PROTON SCATTERING FROM A POINT PROTON 

Although it has been proved [22] that the proton is not a 
pormme particle, 1t 1S convenient initially to solve’ forthe 
cross section for elastic electron-proton scattering from a 
moimt proton. 

The initial assumptions are: 1) the proton has a point 
charge and a point Dirac magnetic moment, 2) the one virtual 
photon exchange adequately describes the scattering process, 


and 3) the electron is a point particle. The Feynman diagram 


describing this process is 
Pe “i 
gq 


Pi Pi 


where Dir Pe and Ps and P, are the incident and final electron 
and proton four-momenta respectively and cfg is the four- 


momentum transferred in the scattering process, i.e., 


fee a= (p,-P,) “= { (p,-p,)*- (B,-E,) *}= (P,-P,)*. 


Brmerder to obtain the cross section for this elastic 
scattering process we make use of S-matrix formulism. The 
S-matrix 1S a matrix connecting the initial and final states 
of the process shown above. The electromagnetic interaction 
responsible for electron scattering may be treated by per- 
turbation theory. The first approximation of the S-matrix 


is given by an integral over the interaction Hamiltonian, 


i 





Ss = ie Vat H(t), 
where H, (t) may be written as a three dimensional iavegral 
over a Hamiltonian density By (erie) 
Voges 
H, (t) = \a a H, (r,t). 


“Since the electrons can be described by a second quantized 
Dirac field and the proton's electromagnetic field may be 
described by a vector potential, the interaction Hamiltonian 


density can then be written as 


H, (x) = - 3) (x)A, (x), 
where x is the four dimensional space time vector. The four- 
current of the electron can be written in terms of the four 
Gemeonent state vector and Dirac matrices as 
J’ (x) = - ie ¥,(x) y,,¥. (x), with ¥(x) = Y"y,. 
pb fe bho 4 


The matrix element of this scattering, Hess is then 
“> 3 
.- =e - |] MY ie a= 
Hey ie \ p OY, , CA, Cx) ie (II-1) 


giving the interaction of the electron current (from its 
initial state to its final state) with the electromagnetic 
field of the proton. By making use of the first Born approxi- 
Mation, the electromagnetic field due to the proton's transi- 
tion current from the initial to the final plane wave state 


can be obtained by solving the inhomogeneous wave equation 


(9% 2) aya = am celal > 
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where the expectation value of the proton's current operator 
meethe transition current creating the potential A (x). The 


Senution is 


Aa 

A, (x) = ee Pel(x)YP, (x), 
q 

where (x) is the state vector of the proton. Substituting 


this result into Equation II-1l the matrix element is 
He = - 4ne"\¥ (x), ¥, (x) 4S Fp 7,9, (x) ax (mee) 
ical f£ oa 2g a a 


Because of the use of initial and final state plane waves, 
the integration over space results in the delta-function 

6 (Dp +P, - piP,), expressing momentum conservation. It is 
interesting to note that the ee term can be considered to 
be with either source current. Thus the proton may be con- 
sidered as being scattered from the electron's field or the 
electron may be considered to be scattered from the field of 
the proton as it was above. 


The cross section may now be obtained by calculating 





GW) 
dao tel 
—— = 2b : (II--3) 
dn. ey flux 
states 


and averaging over the initial spin states, where 


oie = lH, 16(E,-E,) 


and the flux is the incoming flux of electrons. The sum 

over final states is really an integral over all possible 
proton momenta and over that range of electron momenta which 
are within the solid angle of the detector. These integrations 


over momentum space and the averaging over the initial spins 


aa. 





and summing over the final spin states are tedious but then 


completed Equation II-3 yields 


oa 
bee 
eb 
N 
= 
N 
PS . 
Gh 
5 
N 
Ny]@ 


4 2 ‘ 
ae | & e os 9/2 Ji A ) 
dao - Ly epcoseyes =) \ a ee (II-4) 
dn Riicac 4 E¢ sin e/a) _ r 2E,sin O/2/, 
proton mec? 
p 

where ES is the energy of the incident electron and @ is the 
angle through which it is scattered. Both are measured in 


the system where the proton is initially at rest and 


eateca = 4m. E,sin*/2 


Eg Beomeche final electron energy given by 


E: 
1 


£ 2. 


ie - sin°0/2 





mc 


In the limit that mc /E approaches zero the first 
enclosed term on the right of Equation II-4 was calculated by 
Mott [5] in 1929, using the Dirac equation, to be the cross 
section of an electron scattered from a coulomb potential; 
in the case of a proton that physically means that the proton 
is a point charge and that it has zero spin and an infinite 
wees. , it doesn't recoil. Let us cali@ehat the histori— 
Gal Mott cross section. The final enclosed term on the right 
is the recoil term which accounts for both the kinematics of 
the recoil and the fields associated with the accelerating 
proton. When the historical Mott cross section is multiplied 


by this recoil term it is generally called the Mott cross 
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section. That nomenclature is followed in this work. The 
paneo/2 part of the second enclosed term accounts for the 
scattering Cee by the Dirac magnetic moment of the 
proton while the 1 in that term accounts for the point charge 
scattering. Equation II-4 thus provides the cross section 
foreen electron scattering from a point proton with a point 
Dirac magnetic moment. However, there are two corrections 
which must be made to Equation II-4 in order that the theo- 
rectically predicted cross section of a proton agrees with 
Memisucenents. 

The first thing which is known to be wrong with Equation 
II-4 is that the magnetic moment of the proton is not equal 
to one nuclear magneton as is predicted by the Dirac equation, 
but is measured to be Bis = + 2.79275 + 0.0003 nuclear mag- 
netons. Also the prediction of the Dirac equation of zero 
magnetic moment for the neutron compares badly with the 
measured neutron magnetic moment of bn =)=159135 510.0002 
nuclear magnetons. 

Pauli [23] has shown that the Dirac equation for a 
Gm@emgead particle interacting with an electromagnetic field 
can be modified so as to represent a particle with an arbi- 
trary magnetic moment by addding the operator 

= Joye ape } 
v LU 
to the equation. Foldy [24] has verified that such a term 
satifies the restrictions of being Lorentz covariant and gauge 
invariant and has shown that the general interaction terms 


form an infinite series involving arbitrarily high derivitives 
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of the electromagnetic potential evaluated at the position 

of the particles The series of coefficients of these terms 
can be interpreted as a series of moments of the charge and 
current distributiors associated with the particle. We will 
mecurn to this aneemenen later in this section. 


Beeause the four potential, A,, can be regarded as the 


ji 
potential produced by the electron, the Pauli term can be 


shown to reduce to 


K 
—k_o ( _ ) 
Pp. p 
am? i TN ry 


which is added to the vu term of the nuclear current. Herne 
as = 2.79- 1.0 = 1.79 is the anomalous magnetic moment of the 


proton in units of nuclear magnetons and 


oe aes 


= OIE 7 
7 2% 7 vow 


fo) 
bv 

When the Pauli term is inserted into the matrix element of 
Pouation Ii-2, the cross section for a point Dirac particle 


with its normal point magnetic moment and now also its 


anomalous point magnetic moment is obtained: 


2 
22 (a2) i ee ee 2 
a ee {1+ man Z(t) tan (jr 
o* point POEL 4am®c P 2 *p , 
Proton 
where 
4 2 
(22) = ( en eos Se ) 
eee AE? sin*e/2 7 2Esin?0/2 
= ae 5 
mc 
Pp 
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The second change which must be included in the description 
of scattering from a real proton will now be considered. In 
1955 Hofstadter and McAllister [22] measured the cross section 
of the proton and showed that the structure of a proton is 
more complicated than a point’ charge and point magnetic moment; 
the proton scattered as though it had a charge and magnetic 
moment "cloud" spread out to an rms radius of about 0.7 fm. 
Subsequent work (1960-61) on the proton at Stanford [25,11] 
showed that the magnetic moment associated with the Pauli term 
had a softer or more spread out distribution with an rms radius 
emeabout 1 fm, while the distribution of the charge had an rms 
radius of about 0.8 fm. These rms radii were calculated on 
the basis of Equation I-3 by assigning proton form factors 
Fy (a) to the Dirac charge and magnetic moment distribution 
and Fo, (a°) to the Pauli magnetic moment distribution. The 
addition of the, anomalous magnetic moment to the interaction 
and now the provision permitting the Dirac and Pauli components 
to have spatial extent changes the proton contribution in the 
matrix element of Equation II-2 from 

Pg (Xx) 1%; (x) 


to 


(By Pg) Fa, (a7) )e;- 


7 2 ae 
G 

By BY od UU 

The elastic electron-proton cross section with this 
modified current matrix element has been calculated by 
Rosenbluth [26] to be 
do = (22) 2 —t f(r, Ne 2,0 evel 
aA an ae a ie ree 5 4 2\F p* HF ap tan (5) 78 Pop (II-5) 


and is called the Rosenbluth cross section. 
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B. NEUTRON FORM FACTORS 
The neutron has zero total charge. Both the electron and 
neutron have eameuie moments, however, and thus they can 
interact through the magnetic dipole-dipole interaction. 
Also, an incident electron has a magnetic field associated 
Wee it becausewof its charge and this will interact with 
Magnetic moment of the neutron. These interactions are well 
understood and must be accounted for in the calculation of 
the electron neutron scattering emer 
There remains a third interaction which is generally called 
the neutron-electron interaction in the literature. This in- 
teraction is expected to be both spin and velocity independent 
if there is a separation of electrical charge in the neutron. 
Thus an electron, or any charged particle, penetrating the 
extended charge distribution of the neutron will be subjected 
to electrostatic forces. These may arise from two sources: 
1), the intrinsic charge separation in the neutron analogous 
to the spreading out of charge in the Dirac proton and evi- 
denced by the form factor F 


analogous to F of the proton, 


in ip 
and 2), the anomalous magnetic moment Of “Ehe Wweutron, des- 


@erbea by a form factor’F analogous to F for the proton, 


an 2p 
if the anomalous magnetic moment has a spatial extent. 

Foldy [27,28] has shown that since the neutron satisfies the 
Dirac equation, and since it possesses an anomalous magnetic 
moment, then some charge separation is to be expected. He 
compared the amount of charge separation predicted because of 


the neutron's anomalous magnetic moment with that charge 


separation inferred by the very accurate thermal 
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neutron-clectron measurements of the neutron's attractive 
motential for the electron made by several groups [{16,17,29, 30]. 
The predicted and inferred charge separations agreed within 
Beperimental error, implying that the so-called intrinsic or 
Dirac separation is not present. Since this requirement for 
charge separation within the neutron simply because it has 
an anomalous magnetic moment is not intuitive, this will be 
examined in more detail below, as will the results of the 
experiments that measured the neutron's attractive potential. 
Several methods of measuring the attractive potential 
have been used. Only the first [16] of these is outlined 
here. If one measures the coherent scattering of neutrons : 
by spinless (in order to remove effects due to magnetic in- 
teractions) atoms containing a nucleus and Z electrons, then 
the coherent scattering amplitude for neutrons whose wave- 
length is long compared to the size of the atom will be the 
algebraic sum of the scattering amplitudes due to the nucleus 
and the Z electrons. By decreasing the neutron wavelength 
femal it is comparable to the size of the elieectronic cloudso: 
the atom, destructive interference between the scattering from 
the individual electrons within the cloud of electrons can be 
achieved and the total coherent scattering cross section can 
Be made to approach that of the nucleus alone. By separating 
these scattering effects, one obtains the so-called neutron- 
electron interaction. The result of the thermal neutron- 
electron experiments is expressed as the volume integral of 


the potential of the neutron-electron interaction. Plots of 


Za 





the neutron (and the proton for comparison purposes) charge 
meaicity and the interaction potential are: givenerne 1 aime 


m—1 (31). 





Py NEUTRON 


Po PROTON 





Figure II-1 

Neutron and Proton Charge Densities and Taeees cies 

Potentials Between These Particles and an Electron 

The neutron-electron potential is denoted as Ovi that 
is, it is the difference between what the neutron potential 
would be if it had no charge distribution and what it would 
be if it had a positive center and a negative outer region. 
This concept of a positive center: and a negative outer region 
is justified by the results of the thermal neutron~electron 
scattering experiments and is also consistent with meson 
theory. The volume integral of this charge distribution 
does not permit an independent determination of the depth 
and volume of the potential. It does however permit the de- 
termination of the r? of the charge distribution by the 


rms 


following calculation. 
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Let 
= \r?p( p Gigs 56 


and then use Poisson's equation to obtain 





(II-6) 
rms Atte? 


However, following Fermi, all experiments have reported 
meaeults in the form of a potential Ve which was the average 
6 which the electron sensed. This required that some aver- 
age radius of this potential be adopted. By convention, that 
Yadius is the classical electron radius (2.81794 £m), which 
implies that classically the bound electrons in the atom 
were assumed to be interacting almost uniformly over the 
neutron, because as will be seen, the neutrons largest Geta 
@eeless than 1 fm. That is, the electrons are not just 
sensing the outer charge layer as might be implied by the 
Moot cal picture of a point electron of low momentum pene- 
trating only the outer shell of the neutron charge distri- 
bution. This is perhaps best seen from the virtual photon 
scattering picture. The wavelength of the virtual photon is 
very long compared to the neutron diameter and the electro- 
magnetic fields will be effectively constant at any time over 
the extent of the neutron with the greatest difference in 
these fields occurring between the extreme edges of the 
neutron. Thus, while it is the charge distribution near the 
edges of the neutron which show the deviation of the neutron 
from a particle with zero charge density everywhere, it is the 


average charge density, and through Poisson's equation the 


Zo 





average potential within the neutron, that is measured in 
these thermal neutron experiments. Using this average po- 
tential the volume integral of the interaction potential is 
then written 


ig ne Z 


e e 

3 Z 4 3 
\ war = v.\ Ais sels Zig 3 TT xX (IL-7) 
e) Oo mc 





where re is the classical electron radius. 
Since Equations II-6 and II-7 are simply different ways 
of performing the volume:-integral they must be equal, there- 


fore we have: 





5 2V 
r = aes ) 
rms 
mc 
thus 
2 
1g eee ee 
rms Z 
mc mc 


where ts is put in the conventional units of million electron 


volts. Then 


2 


» Z 
eee 31.078 we (im ~ (II-8) 


The presently accepted value of Los is [4]: 


Mes = 3720 + 90 electron volts . (II-9) 


This value of bs when substituted into Equation II-8 gives 


2 
rc charge of neutron = - (0.34 Sa. the negative sign 


rms 

showing, since the charge on the outer shell is weighted by 
2 : ; : 

© more heavily than the inner shell, that there is a charge 


separation in the neutron, and that the outer regions of the 


neutron are negative. 
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It wasyState@gaweve that the two facts that tie iieutnor 
satisfies the Dirac equation and that it has an anomalous 
magnetic moment omen sutiierent to aeeountycompletely “for jane 
measured value of the potential V[27, 31,32]. Let us now see 
wy that statement is correct. 

The relativistic Hamiltonian for a neutron with the 
anomalous magnetic moment as attributed by Pauli can be 


written as 

H = BM + @P-K (e/2m) [BO-H -_ipe-5] 
where KO is the anomalous (and total) magnetic moment of the 
neutron measured in nuclear magnetons. This can be reduced 
to the corresponding nonrelativistic Hamiltonian by the method 
of Foldy and Wouthuysen [33] which yields 


2 e ee 
H = pM + Bp /2M - x6, ($5) p div E 


n « (-25)65- [P ae =) Boor) 
4M 
with terms of order greater than (1/M) 7 being neglected. 
Meaethe coulomb field of a point electron located at x 


the above Hamiltonian becomes 


H = BM+ (BP“/2M) + ani (e°B/am*) 6 (x-x +... 


[mere the term containing the delta function is in the form 
of the electron-neutron interaction. Equating the volume 


integral of the neutron-electron interaction potential to 


this term one has [27] 


Vo = TH, (e7/m*) [ (40/3) (x2)) 1 = - 3900 ev. 
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This compares very well with the experimental result given in 
Bewation II~9. Since the neutrom charge distribution cane 
completely accounted for by its anomalous magnetic moment this 
means that there is, within experimental errors, no Dirac 
pemeicle intrinsic ene distribution unless relativistic 
Semeributions fortuitously cancel this Dirac intrinsic charge 
cas tribution. 

The distribution of the neutron anomalous magnetic moment 
was measured by Yearian and Hofstadter and Bumiller [34,35] 
in 1958. They used Rosenbluth's proton cross section as given 
imeeguation Di~-5 withemodifieatdonseto obtamm the meutzon 
Gress section. Based principally on the thermal neutron work 
discussed above, they set Fin = © to describe the nevweron s 
lack of intrinsic charge and magnetic moment distribution and 
allowed F to describe its anomalous moment. The ao 
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Gisoscs section was then 


2 
2 Cle ce El Ge Se | 
dn dn ‘i Zee oT Am? 2 


Mot 


Electrons were scattered from a deuterium gas target at suffi- 
ciently high energies and large angles to enhance the magnetic 
scattering from the nucleons. By using the rather good approx- 
imation for inelastic scattering from a deuteron, i.e., 

Ss = 0 + OFF eney were aole to Cbtain tic... GOfethe die 

d p 7) rms 
malous magnetic moment of the neutron and concluded that it 


was between 0.8 and 0.9 Fermi. 
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Ce CHARGE, MAGNETIC AND ISOTOPIC. FORM EACTORS 

Thus far the four form factors used to describe the dis- 
tributions of charge and magnetic moment of the Dirac particle 
and the Pauli moment have been introduced in the historical 
manner in which they evolved. They may be summarized as 
follows: 

ite is the form factor associated with the distribution 
i charge of the Dirac proton and also with the distribution 
Semchne intrinsic (Dirac) magnetic moment. In the static 
fo = 0) limit Fy, (0) = 1. 

ss as the fornmetactor associated with the anomalous 
(Pauli) magnetic moment distribution of the proton. In the 
eetac limit Fo, (0) = Shs = Oe 

EA iS the form factor associated with the distribution of 
Smeege Of the Dirac neutron and also with the distribution of 
the intrinsic (Dirac) magnetic moment of the neutron. In the 
Seaeic Limit Fan 60) = Oe 

Fon is the form factor associated with the anomalous 
(Pauli) magnetic moment distribution of the neutron. In the 
Stetic limit F,,(0) = © = 1.91. 

It has been shown above that the anomalous (Pauli) mag- 
netic moment of the neutron causes some charge distribution 
Shin the neutron. Therefore the form factor Fou nou only 
tells what the anomalous magnetic moment distribution is, but 
wears Something about the charge distribution. 

The work of Yennie, et al. (1957 [31], Ernst et al. (1960) 


[8] and Hand et al. (1963)[19] introduced linear combinations 


a 





me the four Dirac and Pauli form factors for both Eheoreticall 


and experimental reasons. These combinations were the 


following: 
Z 
G_={F,-—s,uF 
aN lp heer Pp 2p 
= oC 
“m,, ys + p' 2p 


= K 
Gn ale cs one 


In Ref. 19 it is pointed out that these linear combinations 


Goerespond to zero, G and one, Ga unit of angular momentum 


sed 
transferred along the direction of the virtual photon ex- 

changed in the scattering process. Reference [8] added more 
physical content to these linear combinations of Fy and F. by 


showing that the rms radii of and charge and magnetic moment 


spatial distribution zee derived from Equation L=-3 by 
es «| and ls e |" 


at q? = 0. For this reason the descriptive names for Gr and 


Gy are the "charge" and "magnetic" form factors. Additionally, 


Ref. 19 points out that the cross terms of the charge and mag- 
mecic form factors do not interfere with each other in the 


Pesenbluth cross section when it is written in terms of Gp 


and Gu while the Dirac and Pauli terms are mixed into the 


Rosenbluth cross section as products of each other, as can 


been seen by examining Formula II-5. Written in terms of Gr 


and Guy the Rosenbluth cross section for the proton becomes 
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2 
Sc 2 
ee (co) | Ee ee ee ee oe oe 
| we ee: 2 M 
. 1+ + po —— P 
2 2 
4m 4m 
p Pp 


with an analogous equation for the neutron cross section. 


Gp and Gu are obtained experimentally by measuring the cross 


section atat least two different scattering angles, at the same 


a . from the intercept of a plot of G7 versus the coefficient 


oft Gi one obtains Gh and from the slope one obtains Cu : 
p p p 
The static ee = 0) values of these charge and magnetic form 


factors are: 


Sr (gate OE 
P é 
2 
Gu a =O) = 2679 (IT-11) 
p 
2 
Ga (Cae 0). =>0 
n 
CG ne ee 
Bas Phe, 


We now consider what the slope of G tells about the 


E 
n 
charge distribution of the neutron by studying the primary 

scattering event of interest in this work, the e-d scattering. 


This process may be described by the Feynman diagram 


e . d. 
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where eC. and e, are the incident and scattered clectron 
states, d. and d are the target and recoiling deuteron 

states and q is the four-momentum transfer. The deuteron has 
a diameter of about 4 fm; in this work, with electron energies 
of about 35 to 100 MeV, the virtual photons had wavelengths 

in the range of 10 to 30 fm. With these long wavelengths 
Compared to the extent of the deuteron and neutron, only the 
outer charge distributions are measured. It is important to 
note that the thermal neutron experiments [18,4,16,17] give a 
positive slope for 3 


dq 


M0721). This means that the quantity - z ane is positive, 


3 as do two of the e-d experiments 
n 

i.e., that the outer part of the neutron is negative. This 
is so because the square of the rms radius is obtained from 


the expression 
(ee) 
2 2 2 
a = \ meee4inr p(r) dr) 


e 
s 


co 
where the total charge of the neutron = \ Ane aoa a) 
Oo 


thus a negative ne 


ne shows that the charge density, P(r), 


of the neutron must be more negative at larger radii. 
feetorically, theoretical attempts have been made to explain 
nucleon charge distributions in terms of a third set of nucleon 
mene Lactors which are called the isotopic form factors. These 
isotopic form factors are useful because the structure of nu- 
cleons is determined by the strong mt oem lon and caliculations 
can be arranged in such a way that one can take advantage of 
the conservation of isospin. As the form of the nuclear elec- 


ric current operator [36] 
= 1+ T 
Jy (x) = ep(x)¥, z= # (x) 
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shows, the current has mixed transformation properties under 
metcationsine isospin Gpace; mart of it transfomms as a scalar 
ema part as a third component of an DSsevecter. AS Vamconse= 
quence of these mixed transformation properties, the proton 
and neutron current elements may be expressed in terms of an 


meevector and isoscalar current, i.e., 


once) a (x) ue oh x) and 


pes) = poles) eet (x) 

mine eGllectric and magmetic isotopic form factors then 
enter the cross section calcul@tion as the coefficients of 
the scalar and vector current operators in the matrix element. 


The isotopic form factors, where the superscript denotes 


scalar or vector, are: 


ial = 3 
G,=% [Gg +6, | gag % 
p n 
V a al 


Q 
i 


1 — 
E” @ Ls “Ty q*-0 % 
(p= 25) 


Hi 


B= | | 
L 1 ', 
Gry 5 “M,,” eM q?o 0 4 oan 


em eed 
3 
7 [Sy cy | me 


we 


(HoH) 


ieere have been attempts to give the” isotopic form factors 
further meaning. While these attempts have not been completely 
successful, it is useful to consider them briefly. A state of 
mW mesons, produced by a photon, must have zero total charge 
emer be odd under charge conjugation. Based on this, Drell and 


Zachariasen [37] have shown that only an even number of pions 
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Sentribute to the isotopic vector form facton and» that only 
an odd number of pions can contribute to the isotopic scalar 
form factors. They also have shown that charge conjugation 
does not permit a single pion to be produced since the photon 
is odd under charge conjugation and the mo, the only ™ which 
conserves charge, is even under charge conjugation. The 


possible intermediate states then include: 


Wl 
ane 
bo 
ee | 

ate 

: 

| 


ee oe ue 





Figure II-2 


Possible Intermediate Photon-Proton States 


In the diagram in Figure II-2 the possible existence of 
pion-pion forces were disregarded; but we know that there are 
resonances in the two-pion and three-pion systems, and we 
expect that such resonances may have an effect on the elec- 
tromagnetic structure of the nucleons. Chew et al. [38] and 
Nambu [39] have shown that such propagators as depicted in 
Peroure 11-2 (b) and (c) do in fact contribute to the form 
factors in such a way that they have poles at the masses of 
the Po: Ww and ie The rho (isospin = 1, spin = 1) coneriouce> 


to the isovector and the omega and phi (isospin = 0, spin = 1) 
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eontribute to the as0scalarm formetactors. On thesassun btion 
met Virtual clouds of mutually strongly interacting bound 
mrons are responsible for the electromagnetic properties, 
these clouds show up -as resonances in dispersion theoretical 
expressions for the form factors, for instance 
cpu eee (clare 
Gelde) Sq \ aD (I-13) 
m Sig aE 

where m corresponds to the threshold for resonance creation. 
Furthermore, if one assumes that these vector mesons (spin =1) 
appear as sharp peaks in the spectra of weight ene wires icean 


i.e., neglecting the widths (P(100 Mev), (10 Mev), (3 MeV), 


Peation IIi-13 becomes 


S 
sc) = [-5,)* 29] . 
ames msn ee (aes) 
1+ 45 
M 
S 


which is the form factor expression which Clementel and 
Villi [40] have proposed, where oAii is a constant. The term 
1-S.3 represents the possibility of high mass resonances, 
i1.e€., it represents those resonances which have such high 


masses that their poles may be replaced by constants. M,_ is 


S 
the mass of the scalar resonances. There is one term for 
Gach resonance. Similar expressions for Gr, cs and cs were 


likewise developed and were used by deVries [10] to fit proton 
and neutron data. This fit is discussed in Appendix A and 


provided the normalization for the proton data obtained in 


emis work. 
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Earlier in Section II it was observed that the outer 
megions of the neutron have been measured to be negative in 
the thermal neutron work. There are several theoretical 
reasons to believe that the outer regions should be negative 
[32,38]. Among these is. the fact that the anomalous magnetic 
moment is almost entirely a vector in isotopic spin space, 
i1.e., the anomalous moments of the neutron and proton are 
nearly equal in magnitude with opposite signs. Additionally 
Since the isovector meson, the rho, is the least massive of 
the rho, omega, phi mesons, the uncertainty principle as well 
as the Yukawa potential indicate that the outer regions will 
be dominated by the rho meson. 
fe METHOD OF EXTRACTING THE PROTON ELECTRIC FORM FACTOR FROM 

fae PROTON CROSS SECTION , 

The Rosenbluth cross section for the proton (Eq. II-10) 


Sam be written as: 


2 
on 
ee) (a2) [5 Be Dien? | E 
(2 = \dn mo eee ee (au te) 
ioueo FON Mott p p 


where T = q°/4it0 This can be simplified by using the scaling 
law proposed by Lehmann et al. [41] and others [19], which 


states that G 


M 
Go = ee : (II-15a) 
a Pp 


Schumacher has shown that we should not expect isoscalar 
form factors to equal isovector form factors [42], and that a 


set of scaling laws which relate the form factors Guy and Cu 
Pp n 
to Gr and Gp are 
p n Gy = HGR + HSE 
p p n 


AO 





and 


Ci Set HS 
n Pp n 


For our experiments, magnetic scattering generally accounted 
for only a percent or two of the total scattering in most runs 
and at most 10% of the ee Fhig in a few runs; thus the dif- 
ference in these two scaling laws, which was at most a 1.0% 
difference at owesa-, then made a difference in the cross 
section of less than a nominal 0.01% up to at most a 0.05%. 
Thus the conventional scaling law was used to reduce the cal- 
culations. en the magnetic contribution is accounted for 
by using the scaling relationship, Equation II-l5a, Equation 


II-15 reduces to 


p SE, 
oo £AD oa 
= ao) = Tee [1+ gil ot (147) (27) we tan” 5 |, (II-16) 
where p 
0 ie 
an yee SHO) 
Mott 


Because the last two terms in the bracket are due to the 
magnetic form factors it is convenient to group them together 


by defining 


2 Zz 2 6 
Creep ZI (itt) ten = . TI-17 
< ma) ; ) 
Because the proton cross section can be written as 
doP ly 2 
a = Om (OG (II-18) 


where oe is the proton form factor, we can obtain the proton 


Seeetric form factor by combining II-16, IIt-17, and II-18 as 


aa? 
2 (14+T P adn 
Ge = ge (1tT)  _ We? - KP da (II-19) 
14+C M ? 
eS any Pp oF (0) 
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where 


ee eee i ie 
M (1+C,,) 


is defined as the proton magnetic correction term. Thus by 


measuring the proton cross section and dividing it by oF (8) 


P 


M one extracts the 


and then multiplying this quotient by K 


Seemacre Of the proton electric form factor. 


ie LHE DEUTERON GROUND STATE 

Before the elastic electron-deuteron calculation is studied 
it is useful to consider the ground state of the deuteron. 
Maes 1s particularly true in this case of elastic scattering 
because the deuteron is in its ground state both before and 
after the scattering event. It is an experimental fact that 
the deuteron has but one bound state and that is the ground 
state. Since the discovery of the deuteron in 1932 many of 
its characteristics have been measured and provide useful 
guidance in interpreting the electron deuteron scattering data. 
Some of the more important are discussed below. 

1) The neutron and proton are bound together by an energy 
eee 256 + 0.0015 Mev. 


2) The angular momentum is 1’, and the parity is even. 


3) The magnetic dipole moment is i OrSS7 59570, 0000s 


~ 


nuclear magnetons. 
4) The electric quadrupole moment is 0.00282 + 0.00002. 
5) The root-mean-square electromagnetic radius of the 
Beweeron is 1.96 + 0.02F. 
6) The deuteron potential has a repulsive core with an rms 


radius of about O.4F. 
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7) The mass /ofsthe deuterons .187 s75souMey . 

Elementary calculations solve for the deuteron wave 
mametion in a square well neutron proton potential using the 
measured binding energy. The solution shows that the neutron 
and proton spend most of the time outside of the square well 
potential. The ground state is mostly a IS"state;i.e., an 
£= 0 state, which is indicated not only from general energy 
Sensiderations [43], but from the fact that the algebraic sum _ 
of the proton and neutron magnetic moments nearly equals the 
deuteron magnetic moment. This indicates that the intrinsic 
spins of these particles are aligned and that there is no 
Orbital angular momentum. If the nuclear force were spin 
independent, one should expect to observe a singlet state, 

FS, with the same energy. However, no bound state with J = 0 
has been found in the deuteron. This indicates that the 
neutron-proton force is stronger when the spins are parallel 
than when they are antiparallel; however, this does not mean 
Bie the antiparallel or singlet force is zero. In fact, the 
Singlet force is almost strong enough to produce a bound state. 

Although the —e description of the deuteron is accurate 
to first approximation, a better description is possible. 
First, one notes that the eesti sum of the proton and 
neutron magnetic moment is 0.87925 + 0.0004 while the measured 
magnetic moment of the deuteron is 0.857393 + 0.00003 nm. 

This difference, about 2%%, is larger than the experimental 
error. The second problem with the description given above is 
that a pure S state implies that the quadrupole moment is zero 


since the wave function is spherically symmetric, while the 
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Geuteron in fact has a positive quadrupole moment. The only 
other state in which two particles give the same angular mo- 


5) 
mentum and parity is the D. state (£= 2, S$ = 1, J = 1); 


sl 
The conclusion is that the total wave function of the deuteron 
must have the form ¥ = al. a + an where the subscripts indi- 
cate S state and D state and the coefficients satisfy 


ss + at = 1 in Onpdersato monpmalize the total wave function. 


The exact value of any that produces a quadrupole moment in 
agreement with the experimental value depends on the radial 
Mart Of the function ae which in turn depends upon the de- 
tails of the radial dependence of the neutron-proton potential. 


Values of the D state probability, Py = aX, 


in the range of 

3 to 7% have been obtained from calculations using different 
potentials [44]. A mixed deuteron ground state with = AN 
produces an electric quadrupole moment and a magnetic dipole 
moment, both in qualitative agreement with observations. It 
is one of the objects of this work to help determine which 
wave function best fits these values and from this determine 
more properties of the nuclear potential. Different wave 
functions which are considered in this work have different 
percentages of D state. However it is possible to keep the 
Same percentage of D state and to change the quadrupole moment 
by changing the range of the tensor potential. This potential 
cannot be spherical as that produces a pure S state. The 
positive quadrupole moment indicates that the two nucleons 
prefer to line up with their spin vectors along the same line 


rather than side by side. The tensor potential which is a 


function of the angle between the spin vector of the neutron 
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and the proton and the radius vector r separating them accounts 

for this preferential alignment of spins and is a potential 

which must be added to the dominant spherical potential. 
Several textbooks [43,45,46] cover in detail all the scon— 

cepts that were briefly discussed above. Having summarized 

the deuteron ground state we now consider in detail the 


elastic electron deuteron scattering process. 


fe ELASTIC ELECTRON-—DEUTERON SCATTERING 

Jankus [13] in 1956 derived the cross section for elastic 
electron-deuteron scattering in the first Born approximation. 
In his calculation the nucleons were represented by point 
charges and point magnetic moments. The electron was con- 
sidered extremely relativisitic, but the deuteron was treated 
Menrelativistically. 

While the nucleons have spin one half and two form factors, 
the deuteron, with spin 1, has three form factors. This is 
consistent with the work of Glaser and Jaksic [47] and others 
[48] which predicts that the cross section for scattering a 
Belativistic electron from a potential of spin J will contain 
Zoer 1 form factors. 

Fankue's result, written in the more current terms [49], 
eae including form factors for the proton and neutron, i.e., 


~ 


permitting Jankus's point nucleons to have extent, is 


da 
ac. Cs feet. Sa 2eere. 2ng? (q2)[ cae) 
pe (OM (9) (Sg (a J+ Gn Ge la )+ Sue, (aq) | 142 (14m) tan™ > jy, (21-20) 
2 2 2 2 
where Eq /4M 5 and Go @ Va; Gjlq ) and Go(q ) are the charge, 


quadrupole and magnetic moment form factors of the deuteron. 
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These are given by 


Go = (G+ G., )Cps (IT-21) 
p n 
where 
Cae tee \ |w? (x) + w(r) | (g =) dr (II-22) 
E ra oO 2 
mer che charge from factor; 
Go] (ea Go Ce (r= 23) 
2 as Ee A 
where 
2 
= i ‘[p@we - 2) hh Dae (11-24) 
mer che quadrupole form factor, and 
2 { Jeg+( }e,} 
Gy = 2 va je, + Gu Cot Gp + Gp C,., (II-25) 


Pp p nh Pp n 


Beeethe magnetic form factor, where 


s s {2 (x) - $20 ]5, cof +a weds ig E 5)dr (11-26) 


ee 


and 


oe 3 \. w? (x) 5, (a lee cme Jax. (II-27) 


(G,, + G,, )C,. determine the contribution of the intrinsic 


CUM S 
p n 


magnetic moments of the proton and the neutron to the cross 


section and (G.. + Sew es 


of the convection of charge in the deuteron to the cross 


determines the magnetic contribution 


section. 
Here eh and Ma are the mass of the proton sand deuteron 


respectively and u(r) is the S-state wave function and w(r) 
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mewthe D=state wave L£unctronmormalized sco Ehnat 


CRC es w? (x) jax pb (II-28) 
O 


and Iola = and 35 (4 =D are spherical Bessel functions. 


mne static limits of the functions above are 


“a 


Pe 
Gii(O a= Uns 9 CRORES 
M 
seal [ _3 ogg | 


(0) )=Gy (O)=5(H +H), 


| 6, (0)+ G, (0) | —NGere40)= 37 21G, (O)'G ' 
p be S S 


M M 
p n 


2 3 
Cp (0)= 1, Co (0) = M,Q, C, (0)= Jk 5 Py (II-29) 


and 3 
Cy (0) = 3 Py, 
Q is the quadrupole moment of the deuteron, ee and a are 


the proton and neutron magnetic moments, and 
co 
P = | 4 ie Gr 
D 
Oo 


tmomeene D~state probability of the deuteron. - 

Introducing the nucleon form factors to Jankus's results 
of course assumes that the nucleons are not distorted in the 
deuteron. The fact that these nucleons are so weakly bound 
Mac they spend a sizable portion of their time far apart 
makes this assumption of non-distortion plausible. 

Gross [50] and Casper [20] have calculated relativistic 
corrections to the electron deuteron scattering problem. 


While these corrections could come primarily from relativistic 
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modifications of (a) the deuteron wave functions and (b) the 
nucleon current, the latter have been shown to be small by 
Gourdin [71] and also by Gross. The relativistic modifications 
of the deuteron wave function depend not just on the Lorentz- 
contraction factor, which is ectuallveascentraction 1ngive 
directions [31], but also on the fact that the equal-time wave 
function in the rest frame is not the same as the equal-time 
wave function in the moving frame [20]. Gross goes on to 

solve for the corrections to the deuteron wave function and 
then obtains workable eee enone for the corrections in terms 
and G 


of u(r) and w(r). Gross's expression for G are 


ge Sa 2 
then much more complicated than those of Equation JI-21, II-23, 
and II-25. For example, Gross's complete result for the deu- 


meron charge form factor is 


» © ») fo) 
2 2 Zee. ig 2 
Ce aC gq /32m2 )\ (u +w Jig la 5) dr- —<,. ( (u -uru' 
Oo 24M e) 
D) : 
-uku + Iw°- wrw - whw )[59la H+ 35 (a 4) jar 
2 ow 
+ — \ (urg*+ wr" + wiw'@ ) J fa =) ay, (IL-30) 
16M, fe) 


where prime refers to differentiation with respect to r and 


2 
k = r?| ao - a|, a= [Me, and € = deuteron binding energy. 
ae if 


*~ 


dhe first integral is the usual Jankus nonrelativistic result. 

Schumacher and Bethe [44] and Elias et al. [49] note that 
fae relativistic corrections of Casper and Gross [20] to the 
Beweteron form factors, to order Gane 


p! enter simply as a 


EFector (l- Toe on the right hand side of Equation II-20, 
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where T = q°/4M. Equation II-2G6 can be further simplified 
for small a by using the scaling law [42] 


M 


a 
Gy > aa 
p 


(IT-31) 


fech is not a law but an approxXimate fit to the data which 
can be derived theoretically from the nonrelativistic quark 
model and SU6. (SU6 also predicts tain = 1.5 which is not 
quite true.) 

By Making the relativistic correction as noted by 
Schumacher above and using the scaling law, Ecuation II-20 
eecomes 


aot 2 


A See mer 7, Dae a 
ae of (9) {62 (g ) + 9” Go (q ) + Ales G ieee (1+7)) tan 2) Ge. & Cy 


(II-32) 


Rewriting this using Equation II-21 and dropping the quad- 
rupole term as it is less than 5 x yoo at the highest value 


of q? reached in this experiment, we then have 


d 
Ges un) ee 207 {1+ 2 227 (142 (14n) tan” 8) ha 3)%. (II-33) 
Bs 


Meewe define C. = fru? (142 (147 evs ae (I1I-34) 
. d Creal 2 
and note that 
d 
do d 2 
then 
2 2 ws rT 2 
Ga = ee =p Cp (14C 5) (1- >) 2 (II-36) 


Additionally, if we define KC as the deuteron magnetic and 
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melativistic correction we have: 











Z 
Joma i il) _ 
dad l- oY 
me obtain the so-called deuteron electric form factor Gr 
d 
meer Equation Ii-36, i1.¢€., 
Be | ce: 
e2 = (c, +6, )2e2 = ——89 = okt xd 88 | (ere) 
d p n (1+c,) C- 5) 0, (8) 


With both the deuteron and proton electric form factors avail- 


able we now obtain the neutron electric form factor by first 


i 


forming the ratio of G EO- Guess. 4 
d Ep 


Gy (aa) (G, (q5) + 6, (ag) Cp 
a) G. (42) cake” 

Eo dp Ey P 
where the ae and qs indicate that the proton and neutron in 
the deuteron have received a four-momentum transfer different 
from that of the proton which was the nucleus of the hydrogen 
target. This was so in this experiment because the incident 
electron energy was such that the recoil electron from both 
hydrogen and deuterium targets had the same final energy. 

Thus the transmission efficiency of the spectrometer and the 
detector efficiencies were held constant for both targets, but 
the four-momentum transfer was slightly different. Since the 
meoton electric form factor is available from deVries [10]b' 


fit we may scale the Ga (a5) measured from hydrogen to the 
p 


2 ; : : 
Gp (qa) meesuneéd me@deuterium. Thus the electric form factor 


S10) 





mee the proton in the hydrogen that would be predicted from 
ime hydrogen measurement, if the hydrogen had experienced the 
same four-momentum transfer as did the proton in the deuterium 
is given by 

Ga (a4) 


Ze Z eae ot: ¥ 
Gy (q3) = G (q) = : (II-40) 


Z 
P P Sosa deVries 


Poeuetion II-39 may be rewritten in terms of the measured 
quantities, the deuteron and proton cross sections, by using 


Equations II-19 and II-38. This gives 


2 2 a 
OG 2 
Seg)? Sp “Face ft Se oF (9) 
p n : Tr weun. m 

5} = ae \ ae ga 2 (II-41) 
Ee (a) Gets) p cee 
E Pp m K 

p mda 


hew forming this ratio, i.e., Equation II-41, with the 
same four-momentum transfer for the proton in both the 


deuterium and hydrogen we have: 


MMe )+c (q2)\c aot 5 S| i a 
E rs En da E —_——— oe) K E p 
a 


5 = ea BE | 4) PS Gee } 
ee) (a5) do" of p Cn CoD, 
ED 2 dn one) Kn ED a deVries 


which gives 


2 1 PE - 
PE (qq) oF (9) Kau 5 Gh scm) , ; 
G oO 
2h a) a eee deVries 12 


where R is the experimental ratio of the deuteron cross section 


EO the proton cross section. 
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Note that oP (9) and ce were calculated at the incident 
energy of the electron on the hydrogen, which was higher 
than the energy of ine electron on the deuteron, at which 
o° (8) and KO were calculated. 


The neutron electric form factor, G Coy then is 


E 
n 
given by 
2 
d\ *%/(G, (g7) 
oP (9) K ae 
G eo) = R mr mee ee +t _i.ole (Ga) 
14 d D C da 
En oo (ojmEe G, (q<) ee ED 
m m De a deVries 
(=) 
which may be rewritten as 
G 
ao.) = tao Ne a GS | (IT-44b) 
eq’ ~ |G. ¢c ; BE ‘4g’: 
n ED i p 


ce (a4) is the Gp which includes Gross's relativistic 
n n 


@emrections. it appears at first to be quite different from 
those of Stewart and Mader because they used the nonrelati- 


metic Geuteron wave function to calculate Gr and then added 
n 


the change in G which®is due to this welatiwastic correction 


1 
n 


in the deuteron wave function to the nonrelativistically cal- 


culated Gy . They used Gross's equation 10 [20] in a form 
n 


which is:practically model independent within the range of oe 


Gomsidered here [2], i.e., 


2 2 
= A = 
(G_ ) G_ + 4G Gy + ee ‘ 


Ey rel n En n 
Mires can be shown to be very close to the correction of 
Schumacher used in this work by considering the following. 


Seenumacher's relativistic correction enters as a factor 


(1 - m7 Dy onthe rigue hawd side Of Equation 12-20. ‘Thas 
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Soerreétion is carried forward to Equation II-43. If we 


form the nonrelativistic ratio of Gp to Gn we have 


n p 
@ Pp a G_ ( os 
E ohio) eiemale jean “5 
n m in ee 
a = {| R a a 5 —— eS (11-45) 
E non O° (6) KP Gh (qy) E 
Rel a” a 
Define 
a a i 


nC 
iteg KC were inserted in Equation II-44 in place of K, we 


would have 














E 
—n 
G e 
ES Rew 
Since 
G a; (G., (q2) 
Eo oak ae Bo ae 
oy pil of (9) xP Grew =) Cre 1/2) es oe 
pm yRel m ™m gE ‘4g 
Pp 
then: 
2 
G G G Pp a\3(G_ (q_) 
4 (ey mJ 
Mee {in| _{ a) « won ee (Pe 
G, | |G G _ a p 2 2G. 
0 < 
a ED Rel ED non e) Kn oh (qq) = 
Rel 
(II-47) 
and 
7 2 2 
oP (0) KC 2) ore oa) Ga a! 2 
NG — R---_ —_ — ’ (11-48) 
a of (g) xP G. (q2) CE gm- 
m m E d Pp 


where T = a” 4m, has been written explicitly. Since the 
coefficient of q° /8m,, in Equation II-48 is experimentally 
within 2 percent of unity throughout this experiment, then we 
see that Schumacher's method and Stewart's and Mader's method 


are in good agreement. 
as 





GS. JHE ROSENBLUTH PLOT 

Rosenbluth derived the cross section for electron-proton 
scattering based on the one photon exchange. Jankus made the 
Seme assumption in calculating the deuteron cross section. 
Examining Jankus's result, given in Equation II-20, and, re- 


Peoauced below, 


d d 
do do oe DO 20 





(II-48a) 
which can be written as 
a a 
oO Co 2 
= a fa) Ga : | - (II-48b) 
; Mott 
memean be seen that i£f the deuteron cross section is measured 


at different angles but at a constant an menen a plot of en 


versus mean) tan@ 3) should give a straight line with 

Ge + Bn*cs as the intercept and net as the slope. Thus a 
linear plot of such data would indicate that the one photon 
process accurately describes the phenomena. Such a plot has 
been made in Section VB. 

Equation II-48b cannot be used directly because that 
equation is presumed to give an absolute cross section. The 
absolute deuteron cross section is obtained through knowledge 
of the hydrogen cross section. Since deVries b' fit is the 


Meot Lit to the proton absolute cross section, this is used to 


Mermalize the data. That is 


doP _ f aoP gP 
& = (32) = Nae AE) 
Absolute deVries LINAC 


where N is the normalization factor for a given point. 
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Therefore 


p p 
oe ae) . ge) ; (r= 50) 
deVries NPS LINAC 


Because the deuteron and proton cross sections were measured 
under the same experimental conditions, the absolute deuteron 


Mress section is 


aos _ doe 
an nN adn ° 
Absolute NPS LINAC 


He L£Ollows that 


ue Absolute oe deVries uk i LINAC 


In terms of the experimental ratio R, where 


ce (ac?) Age) 
da] wps LINAC “ } NPS LINAC 


aot aaP 
aa mt | i's R. 
Absolute deVries 


aoe aot 2 
an ca igs Gar 
Absolute Mott 


therefore 


Pp da 

(e} 

ae = 2) ae) Re Casi 
deVries Mott 


doP ao 
where (3) ee and ($5) ove were calculated by programs 


TOPROTON and TOPDEUTR which are listed in Appendix E and the 


Bat 


Bacio R,, was the weighted mean given by Equation IV-20. 


Do 





The statvobicel erron 1a Ga was calculated from 


((2 a K G g * aa 
d 
- te Siegen ae lesa Rn 


Ce 


where Sp is calculated stromerquation lV —21- 
m 


Several points were taken in order that the data could 
be checked for consistency on a Rosenbluth plot. The results 


are given in Section V-B. 
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III. EXPERIMENTAL EQUIPMENT AND PROCEDURES 


A. THE NAVAL POSTGRADUATE SCHOOL LINEAR ACCELERATOR (LINAC) 
The NPS LINAC is a 100 MeV accelerator technically similar 
fo the Stanford University Mark III linear accelerator which 
is described by Chodorow etre [51]. A detailed description 
of the NPS LINAC is given by Barnett and Cunneen [52] and a 
Pewematic floor plan is given in Figure IIi-l. The accelerator 
receives up to 22 megawatts of peak power from each of the 
three klystrons. The machine is pulsed at a rate of 60 pulses/ 
fee, Cach pulse being about 1 microsecond in duration. 
Electrons are injected at 80 kev and can be accelerated to 
energies between 15 and 100 MeV with an average beam current 
of about 5 microamperes. The beam energy is selected by the 
beam deflection magnets. The magnetic field of these magnets 
is measured by an NMR probe. The beam is focused down to a 
diameter of about 2mm at the target by the quadrupole magnets. 
The spread in the beam energy is controlled by the slit 


settings. 


B. THE SPECTROMETER AND COUNTING SYSTEM 

The spectrometer is a 10>, 16 inch double-focusing 
magnetic spectrometer described in detail by Oberdier [44]. 
Men detectors span a portion of the focal plane of the spec- 
trometer, thus giving 10 counting channels for each spectro- 
meter setting. Each detector (a NE102 plastic scintillator) 


is 7/16 inch wide and 1/16 inch thick. The 7/16 inch width 
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allows a momentum spread 4p/p of about 0.3% for each counter. 
Behind these ten detectors are two large detectors, both of 
which span the same me Auion of the focal plane of the spec- 
trometer as do the 10 smaller scintillators. An electron 
passing through one of the ten scintillators produces a light 
pulse, which is collected by the photomultiplier tube, ampli- 
fied and transmitted to the control room where it is discrimi- 
nated, and is then sent to the scaler if both backing counters 
simultaneously registered outputs. This triple coincidence 
pulse to the scaler must also occur during the time that the 
master gating circuit is pulsing the scaler in order that the 
count will register. The gating circuit prevents background 
counts from registering between beam pulses. A block diagram 
Meche counting system is given in Figure III-2. 

The scalers in channels 4 through 7 are capable of re- 
cording periodic pulses up to about 150 M, while the other 
Six channels of the ten channel ladder We Plo Teer to scalers 
that can record periodic pulses up to 20 M, - Only data in 
channels 4 through 7 were analyzed since ee. could accumulate 
Meeeewmore rapidly. The count rate in these channels was held 
alee, GO. counts/second one all suns. Countieate cornection 
data were obtained and were used to provide the count rate 
Berrection factor given in Equation (IV-ll)» The correction 
Paetor was linear out to about 100 counts/second, thus main- 
Bening the counts below 60 counts/second was more than 


adequate. 
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The gas target is covered in detail in Appendix D. Both 
Hy and D5 were maintained at liquid nitrogen temperature and at 
eeeut 150 psig in order to have a high density of target 


mucieli. 


fe METHOD OF OBTAINING DATA 

Consider first the data which were obtained during a 
"run" to determine the ratio 0% /oP at a given q? and scat- 
tering angle. The same data were taken for both deuterium 
Pmaehyarogen and consisted of that which was recorded auto-— 
matically and a lesser amount that was recorded manually. 

That which was recorded automatically consisted of any- 
where from 30 to 60 individual units of data, depending on 
the run. Each unit of data was called a "point" and consisted 
of the data accumulated at a given spectrometer setting while 
a given number of incident electrons passed through the gas 
target. In each run, the spectrometer was initially set such 
that the scattered elastic spectrum was well above channel 4. 
The first of the data points was recorded there. Then the 
energy of the spectrometer was lowered between 0.02 to 0.05 
MeV, depending on the energy range the spectrometer was 
sweeping, and another data point was recorded. This process 
continued until the spectrometer setting was below the peak 
of channel 7. The last data point was oo Beeconued., 

After the selected amount of beam charge had been accumu- 
lated at each data point, the following data were automatically 
fed to a teletype machine that both printed a hard copy and 


punched a paper tape, a) the number of counts recorded in 


ol 





each of the 10 channels, b) the number of counts in coinci- 
dence in the two backing counters, c) the time required to 
obtain the given charge, ad) the incident energy of the 
electrons, e) the capacitance of the capacitor used to in- 
tegrate the incident beam, f) the setting of the spectrometer 
in MeV, g) the maximum voltage permitted on the capacitor, 
and h) a number which, when divided by 30 (the full scale 
reading of the integrator), gave the fractional value of the 
maximum voltage at which the integration of the beam current 
stopped. 

The data which were recorded manually included Secondary 
Emission Monitor (SEM) efficiencies and the gas pressure. 

The SEM efficiencies were taken just prior to the first data 
poimt, at the peak of each channel and after the last data 
point. The gas pressure was recorded at the start and finish 
Sie eden liquid nitrogen "fill" (as defined in Appendix C), 
which generally occurred about once at each data point. The 
leakage current in the integration circuits was also recorded 
Manually and was used to determine the ground loop errors, but 
mes emOt used to calculate the cross section as it was always 
fier a Small correction. 

The above data were always taken under conditions that 
were as identical as possible for both the deuteron and 
proton. That is, the experimental geometries, the accele- 
rator beam current, the pressure of the target gases, the 
ecattered electron energy, the detection and electronics 
Systems and all other possible variables were maintained, as 


nearly as possible, the same. 
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Runs were made at values of q? that went from 0.025 fm 7 
to 0.40 fine However, scattering from the target windows or 
a contaminant in the gas masked the data at 0.025 a. while 
high background levels and small cross sections made the data 
taken at 0.40 em" ee Thus the reproducible data 
meoorted here covered a region from 0.05 to 0.35 fm“, 

Senin this range of q*, forward angle scattering was prefer- 
entially selected to emphasize the electric interaction. 

The most significant change in running conditions between 
the deuteron and proton was the change in beam energy made 
so that the scattered electrons would have the same energy in 
both cases. This insured that the transmission of Sac: 
through the spectrometer and the efficiencies of the detectors 
and Ho “data. This @aese 


) 2 
meant that the ne for hydrogen was slightly different than 


would be the same while taking D 


for deuterium. This difference was accounted for in Equation 
II-42. 

The spectrometer was saturated by increasing the current 
meee coils to a setting corresponding to 100 MeV before data 
were taken in either He OE Dd. in order to insure that the same 
hysteresis loop could be obtained on both spectrums. The 
energy defining slits were generally set to give a full width 
at half maximum to beam energy ratio of about 0.4%. The count 
rate was held below 60 counts per second during all runs. The 
number of counts recorded varied from as much as 170,000 to as 


few as 10,000 in a single channel. The number of counts ob- 


tained depended primarily upon the size of the cross section. 
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An empty Larger measurement was made vee tiewleci a co 
trometer setting as well as at several others throughout the 
Seergy range in which data had been taken. Usually those four 
er five readings indicated a "flat" background. Only at the 
lowest ae OfS0R05 moo did any structure show itself in these 
empty target background measurements. Thus with all but the 
“ie =—20205 nme data runs, the full target background measure- 
ments taken at energies well above the peak of each channel 
were averaged and then projected at those same values through- 
out each channel's elastic peak. The q” = 0.05 po? data were 
handled differently because there was scattering from con- 
taminants in the gases and the windows which complicated the 
background spectrum. In that case empty target background was 
mekecneat each and every spectrometer setting of both the hydro- 
gen and deuterium. These cmpty target points were then scaled, 
Meant by point, to the full target background in order to ob- 
tain the full target background under each elastic peak. 

A significant improvement in procedure over other work 
which had been done at this facility is the procedure that 
steps each detector completely across each elastic peak. 
Formerly, the data from all ten channels were used to construct 
one elastic peak by measuring the dispersion between detectors 
and —_ efficiency of all others with respect to a chosen re- 
merence detector. This calculation of the efficiencies intro- 
duced errors around 0.5% into the former method which were 


avoided here. 
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IV. DATA ANALYSIS 


A. GENERAL METHOD OF DATA REDUCTION 

The treatment of the experimental data (see Section III-D) 
is described here. The first step was to change the form of 
the data from punched paper tape to computer card data. This 
was accomplished by a program that required the data in paper 
tape form to be read into an on-line teletype which was con- 
nected to the NPS IBM 360. This program caused the data to 
be punched out on computer cards. Because of a peculiarity 
in the data readout system of the accelerator, decimal points 
were punched as colons. A small program was written which took 
the punched deck mentioned above as input data and replaced all 
colons with decimal points in a second punched deck. This deck 
was then interpreted and checked with the printed teletype out- 
put. Hydrogen and deuterium data were separated and were used 
as input to FORTRAN programs TOP#FORM where # ran from channel 
4 through 7. The TOP#FORM program is listed in Appendix D. 
TOPHFORM was written primarily to calculate the cross section 
measured by a particular channel, but it was also used to 
provide count rate corrected @ta and to normalize all points 
to a standard incident charge, Since the automatic termination 
of the integration at each spectrometer setting did not always 
stop at the same value. TOP#FORM then prints out these count 
mate corrected, charge normalized counts for each spectrometer 
Setting. This part of TOP#FORM was used to plot these so- 


aurrcad Corrected counts, For that plotting run other input 
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data which are necessary to calculate a cross section, such 
as radiation corrections, presSure, etc. were given dummy 
values. These count rate corrected and normalized counts for 
pie rour detectors for both the deuterium and the hydrogen 
were then plotted and the position of each peak and the shape 
of the scattered spectrum was determined by visually fitting 
a smooth curve through the points. A typical plot of such 
data is given in Figures IV-la and IV-lb. The difference in 
energy between each peak and the final spectrometer energy at 
which data were taken was then measured because the radiation 
corrections depend upon AE, and this value was put into the 
FORTRAN program TOPRADCR which is listed in Appendix E. 
Additionally, the peak energy (which is used as the final 
energy of the electron) as well as the pressure of the hydro- 
gen and deuterium are entered as input data to TOPRADCR. 
TOPRADCR has the virial coefficients of hydrogen and deuterium 
and the composition and thicknesses of the gas target windows 
P'Orea aS constants. These constants and the run data men- 
tioned above permit TOPRADCR to calculate the density of the 
gas for each channel and the Schwinger, bremsstrahlung and 
Landau Hadlation Cormections. The rationale behind tie radi— 
meron Correction calculations is given in Section IV-@. The 
average background was calculated for each channel by adding 
all counts recorded prior to the start of the rise in each 
channel and dividing by the number of spectrometer settings 


winch recorded these counts. 
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With these calculations done, the Schwinger correction, 
the product of the Bethe and the Landau corrections, the gas 
Mensity, the SEM efficiency and the background counts for 
each channel were typed on computer cards and entered with the 
appropriate hydrogen or deuterium deck as the data for TOPH#FORM 
which then calculated the elastic cross section and the sta- 
tistical error in that cross section. The method of calcu- 
lating the cross section is given in Section IV-B. The method 
See calculating the statistical error is derived in Section IV-D. 
Four cross section ratios (one for each channel) were 
formed from the cross sections, and then a weighted mean ratio 
Wes Calculated. The electric form factor of the neutron, 


Gp 1 was then calculated by using this weighted mean ratio, Rat 
n 


and the appropriate cross sections and magnetic corrections 
and form factors as shown in Equations IV-14 and IV-15. The 
wWalues of the proton electric form factors and the proton and 
deuteron Mott cross sections and magnetic corrections were 
calculated from computer programs called TOPROTON and TOPDEUTR 
which are listed in Appendix E. Values of D. were tromsdata 


in Appendix C. 


fee dt; CALCULATION OF ELASTIC CROSS SECTIONS 


1. The Experimental Cross Section 


The experimental cross sections have been calculated 


from the expression 


Soo OC 
a Neve eb ‘TT 2 - 

ae. ster 
Exp Nin uPe Igy 


1 (y=) 


So) d 





where 
Nog is the number of scattered electrons counted in a 
Peeecicular channel, 
Nay is the number of fa etaene electrons, 
ure is the number of target atoms per em, 
An is the solid angle defined by the spectrometer 


Slats, and 


= K and ne are radiative corrections (see Section 


b 
IV-B) © 
ny 1S given by 
18 ie No 
"= Rcos @ ! (ave) 
where 
p is the density of the target (gm/em*), 
te is the thickness of the target (cm), 
No is Avogadro's number, 
A is the atomic weight of the target gas, 
and 


~ is the angle between the normal to a plane solid 
target and incident beam of electrons. In the gas target 
analysis is set equal to zero as the angle is not de- 
med, i.e€., rotation of the gas target still presents the 
same number of target atoms per square centimeter. 


~ 


N.. is determined by 
ai 


O CV 
N. s*s—_ (IV-3) 
abi = €fonm 
where 
Q = total charge of incident electrons 
@ = electron charge of 1.60210 x 10°? coulombs 
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© = capacitance of antegrator an geared. 


< 
Il 


voltage of integrator amivolic, 
and 

B — efficiency of SEM, 

Equation IV-l1 is appropriately used to calculate 
absolute cross section, but it was used in this experiment 
even though absolute cross sections were not measured. This 
is a correct procedure because, in this experiment, the ratio 
of the cross sections of the deuteron and proton was formed 
and the systematic errors effectively concelled as all data 
were taken under almost the same conditions for both the 
proton and the deuteron. 

ee the total number of scattered electrons, is 
@ameaimed by integration. Consider Figure (IV-2), where the 
number of detected electrons in a particular counting channel 


is plotted against the spectrometer energy. 


Let E. be a spectrometer setting 


Ni tae number of counts recorded in this channel 
at that spectrometer setting. 


en area under scattering peak 


Ne. = (in counts Slew (IV-4) 
where 
OE = the energy width of an individual detector. 
This is because each detector has a finite dimension 


and therefore counts all electrons in the spread AE. Thus 


iil 





Number of 
Electrons 
Counted 


background 


enn mmr lt 





Spectrometer Energy (E) 


Pugure IV-2. Spectrum of Counts in a Single Detector 


mn obtaining Neg we have: 


iE 


i 
aN 
es = \ aE dE where 


1 


“ is the number of (IV-5) 


counts per interval of energy. This integral is numerically 
approximated by the sum: 
ee ne 
sc I 
The energy width 5E is obtained in terms of E, by using the 


dispersion relation for the spectrometer optics, 


EB De ed 
O 

where D has been experimentally determined to be 3.92 [53], 

Ly is the mean radius of the spectrometer, 16 inches, and 


‘'w' is the width of each detector (NE102 plastic scintilla- 


tors) which is 7/32 inches. 


a2 





Therefore 


& "i 
ee er ee, 


By yaoi ren) . (IV-8) 


Putting these together, we have 


n-L 


N = 2% “aN. +N 


Sc i=1 (E;-E 


a p41) 70 - 0034877 be. 

The energy E. in the denominator is the spectrometer energy 
setting. This setting is very near the center energy of the 
band of electron energies which the spectrometer detects when 
fee tO the reading E.- However, the energy of the electrons 
which hit counter 10 will be higher than any other detected 
electrons and the energy of the electrons which hit channel 1 
Will be the lowest in energy. In order to take into account 
the fact that each detector is recording electrons with dif- 
ferent energies than the spectrometer setting Ess several runs 
were made to determine both the absolute calibration of the 
Spectrometer and the relative energy spacing of each of the 
Bomar Center channels. By comparing the calculated elastic 
peak energy with the spectrometer setting Boe it was deter- 
mined that the peak of channel number 5 was 0.17% higher than 
Phe chergy setting of the spectrometer. Thus channel number 

5 was chosen as the reference energy channel. By successively 
lowering the spectrometer setting from an energy above the 
elastic peak of channel number 4 through sihamuell mnohboyeie i clots 
dispersion of the magnet was measured. This procedure was 
done at LINAC energies of 17 MeV, 30 MeV, 60 MeV, and 90 MeV. 


The average dispersion of these measurements showed that the 


average energy of the electrons which counter 7 (where 7 is 
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4,5,6, or 7) detected was given by 


Ey = E, (1.0 + 0.0039 (7-5)) (IV-9) 


when the spectrometer is set at Ea: 


Finally, 


n-1 
We eh = " 
(24 itMisn [Fi Bj 43) /0-00348778, (1.040.039 (n 5)) CS 


=2 
Ss ) Je ERE oe )(a.831 x 107° steradians) 


€fomM aS 





(IV-10) 
where the solid angle, An is 1.831 x ne steradians. 
2. Corrections to the Observed Counts 
a. Equation IV-9 would correctly measure an absolute 
cross section if: 
iN) Naa the number of counts recorded on the 
scalars, was the actual number of electrons scattered, i.e., 
memene electronic system did not have a finite resolving 
time. Even though the scalars are capable of counting peri- 
odic pulses up to rates of 150 MHz, they lose counts when 
random (Poisson) pulses are counted at high rates. At 60 
counts per second, which was the highest count rate which was 
permitted during this experiment, about 6% of the counts were 
lost. 
asa) N. consisted only of counts due to the 
scattering from the target gas and was not contaminated by 
background counts. 


iii) The efficiency of the detectors were known. 
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b. The count rate correction was determined experi- 
mentally by scattering from an aluminum target at count rates 
up to 500 counts/sec. The number of counts detected was 
pultted as a function of the count rate. Up to about 100 
counts/second the number of counts recorded decreased linearly 


with rate. A linear fit was made to the data in this region, 


and gave _ N(recorded) 2 
Nieount rate cCorrectea, — 1-0.00097N/t (IV-11) 


where N/t is the average count rate determined by dividing the 
recorded number of counts by the time of the integration. 

c. The background was determined by counts obtained 
at energies greater than the elastic peak of either H5 or Do, 
whichever was the target at that time, for each of the four 
Mme-erorc, The average number of counts for each channel, in 
this region, was then determined and was used as the background. 
During each run an empty target background was also taken to 
insure that there was no rise or fall in this extrapolated 
average background under the elastic peak caused by some source 
other than the H, or D.- | 

d. The efficiency of each detector was not deter- 
mined because the object of this experiment was to obtain the 
ratio of the deuteron cross section to the hydrogen cross 
section. Because this ratio was determined for each of the 
four detectors separately, the efficiency of each detector can- 
celled and therefore is not included in the determination of 
the cross section. 

e. An electrometer was used to measure the inte- 
grated incident beam current for each data point. The cir- 


cuitry associated with the electrometer was set to stop the 


TS 


integration and the counting system whenever a predetermined 
voltage across a capacitor charged by the current from a SEM 
was reached. The electrometer did not always stop at exactly 
the same voltage but varied sometimes by as much as a percent 
or two. Therefore both the hydrogen and the deuterium counts 


were normalized to a standard voltage (V ) for each 


Reference 


magi. fnus, if Vea) was the actual voltage at the end oer 


the integration, then the number of counts recorded was cor- 


rected by 
Veeference 


Actual Vaeeaad 


N(corrected for voltage) = N -@av—1.2) 


TMiesGLress section, after correcting Boer count race, 


background and voltage cutoff is 











Vee Ee 
aE “1 i ie Bra) + Nisn Yact BKG 
“ ~ = = = 
ao =i al 0.00097Rate; } (all 0.00097Rate, _,) 
dn —2 
(au ) (eBEeNosems ) (1.831 x 107° specaatenen 
ec A 
SEM 
(Ee -E. _)K one 
{9 ZEA ae ), yee 


0.0034877 (1.0+0039 (1-5) )E; 


C. RADIATION CORRECTION CALCULATIONS 

The proton cross section was derived by Rosenbluth [26] 
with the assumption that a single virtual photon is exchanged 
in the scattering process as shown in Figure IV-3. 

The deuteron cross section, as derived by Jankus [13] in 
the impulse approximation, also assumes that one virtual pho- 


ton is exchanged as shown in Figure IV-4. 
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Feynman Diagram of Rosenbluth Calculation 


Figure IV-3 
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Feynman Diagram of Impulse Approximation 





Figure IV-4 


Unfortunately, when one does an electron scattering 
experiment these processes are not the only ones which occur, 
and the experimenter must take into account competing pro- 
cesses in order to be able to compare experiment and theory. 
In fact the experimental scattered energy spectrum which is 
measured contains the effects of three such competing pro- 
cesses which must be accounted for before one can obtain the 


cross sections for processes as shown in Figures IV-3 or IV-4. 


iy 


These competing processes are: 

1) Photon emission and absorption while the electron is 
gm the field of the scattering nucleus of interest.. The .cor- 
rection to the cross section due to this process was first 
calculated by Seiwineee ean and 1S commonly called the 
Penwinger Correction. 

2) Bremsstrahlung emission in the target before or after 
the nuclear scattering of interest. 

By Excitation and ionization of materzral in the target. 

Each of the above processes causes the energy of a 
scattered delectron to be less than it would be if the only 
process which occurred was the one photon exchange. Thus, a 
radiative tail appears on the low energy side of the elastic 


peak, as shown in Figure IV-5. 


qo 
dadkE 









Radiative Tail 


Ervergy Of Scattered Electrons ~ 


Elastic Peak with Radiation Tail 


Figure IV-5 
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To obtain a cross section for the nuclear se2ctemiug 
event of interest, the experimenter would have to continue 
Gown the radiative tail to zero energy in order to Capture 
Biose scattered electrons which have some probability of 
losing all their energy in the three processes listed above. 
In practice, the measurement is stopped at AE below the 
elastic peak (see Figure IV-5). The value of AE is generally 
Of the order of at least four half-widths (!/2)[55]. Thus 
fewer electrons are detected and the experimental cross 
Mection is less than it should be. Since each of the three 
processes mentioned above contributes to the radiative tail, 
corrections must be made to account for the undetected elec- 
trons with energy less than E, - AE. 

im ochwinger Correction 

mee largest of the corrections is that which accounts 
for photon emission and absorption while the electron is in 
the field of that nucleus which causes the large angle scat- 
tering which is of interest. This correction to the cross 
section was generally about 15% in this experiment. 

The Schwinger correction accounts for many inter- 
actions which May occur in addition to the one virtual photon 
exchange shown in Figures IV-3 and IV-4. Several of these 
Poeeible interactions are shown in Figure IV-6. 

In processes (a) and (b) real photons are emitted 
before and after the nuclear scattering. Processes (c), (d), 
and (e) are called vertex renormalization terms and (f) is the 
vacuum polarization term. Improved radiative corrections to 


electron-proton scattering have been calculated by Tsai [56] 


fee) 








Feynman Diagrams of Schwinger Correction Interactions 


Figure IV-6 


and by Meister and Yennie [57] and include the effects of 
proton recoil. Numerical evaluation of the Schwinger cor- 
rections in this work used the results given in the review 


eecicle by Maximon [58]. — 


where E. and E, are the electron energies prior to and after 
scattering, @ is the fine structure constant, m is the elec- 


tron mass, ae is the four-momentum transfer, and 


A 





x 
—— log (1-t) 


as the Spence function. 
2. Bremsstrahlung .Correction 

When passing through the target material, the elec- 
trons also radiate away energy ence een radiative collisions 
with atomic electrons and nuclei other than the target nuc- 
leus of interest. The energy loss due to bremsstrahlung has 
been studied by Bethe and Heitler [59], who calculated the 
eross section Ee) and the average energy loss due to this 
process. By approximating this cross section Heitler [60] and 


Bethe and Ashkin [61] obtain: 


where Xo is the radiation length, which is the distance over 
which the electron has its energy reduced by a factor l/e. 
In this experiment, the expression given by Butcher and 
Messel [62] has been used to calculate the radiation length 


Xo? 


Ft oan 


40N rz (240.8) p tn(1s3z 1/3) 


where No is Avogadro's number and re is the classical elec- 
tron radius, P is the density in gm/em™ and A is in grams/mole. 
The target, which consisted of the hydrogen or deu- 
terium gas, was in a pressure vessel which had 0.001 inch steel 
windows for the incident and scattered beam. For the brems- 
strahlung correction, and also the Landau corrections, the 


target was treated as a homogeneous molecular compound of 
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Sircel and exzther Ha. or D 


2 ali 
Ho 6: Or Fe, Doo: where the “molecular formula" was calcu- 


5 which consisted roughly of Fe 


feted for each gas pressure. This formula accounts for the 
1 mil windows and the 2 inches of target gas at 150 psig and 
pa K. Following Butcher and Messel the radiation length for 


Suen a compound target was calculated from 





ie 1 ened 
a oe . . . j % 
eo i ee = x" 


In the actual FORTRAN programs the Xo were multiplied 
by the average density of this compound to give Xo at siya 

Bethe and Heitler [59] also give the probability senee 
an electron traversing t radiation lengths will have an energy 
loss less than @® times its initial energy (aside from the 
ionization loss and the scattering from the nucleus of 


Meeerest). That probability is 


y(t /en 2) 
a) ~ T (i+ tf 2) ° 


a 
Hofstadter has integrated this and calculated the cor- 
rection factor to the observed cross section; this brems- 
strahlung correction factor may be written as [6] 


cae ie | Ee 
ae xr tog (ig) * 


~ 


where T is the thickness of the target, including the two 
steel 0.001 inch windows as well as the gas, in gm/em-, and 
E is the geometric mean electron energy while going through 


the target. Ky was about 1.02 in this experiment. 
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So Leonazataon Coimeect ion 

As the electrons in the beam pass through the target 
they also lose energy by excitation and 1oni zation aneene 4.4 
and in the gas target windows. Landau [63] first calculated 
both the average energy a monoenergetic beam would lose, and 
the spread, or straggling, of the energy distribution of the 
Imeem after passing through a thin target. From this distri- 
Petion it is possible to calculate the fraction of electrons 
which are not detected because data is taken only down to 
E,~4E. Figure IV-7 shows a theoretical spectrum of monoener- 
getic incident electrons after passing through a layer of 
material. FV is the full-width at half-maximum. ase 12 tthe 
most probable energy loss due to ionization and EINC is the 


incident energy, and N(E) is the number of electrons detected 


Peewenergy FE. 


N (E) 
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Biftect of Ionization on Spectrum of 
Monoenergetic Incident Beam 


Figure 1V-7 
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Following Landau's [63] method, but using the general 
moetation of Breuer [64], one can calculate the freetion, (AE), 
of the beam one does not detect (because of ionization losses) 
if data taking stops AE below the peak in Figure IV-7 as 
ub ane: 


2 
SGA) \ N(E)dE / \ 


Oo 


Fee 
N (E) di 
O 


Landau [63] gave a graph of ¥(AE) but he also derived 
an approximation for ¥(AE), as have others [65], which has 
been used in this work: 


¥(AE) = x, 


where the value of Was provided by Bérsch-Supan [66] is 


Yn A+A ; 
A+1 


"w= A.(1— 
memory = Y —~ i = -0.4228..., 
Y is Eulers constant, 
X = iB - 0.05 (is an approximation to the A as defined 
by Landau and Bergstrom [67], 


a NZ, 
©. 54 = n,A, i MeV, and 


um 
I 


T = target thickness is aoa - 
n. is the number of atoms of element (i) per molecule 
of the target which here was considered to be a homogenous 
mixture of the 0.001 steel windows and the gas in the target 
vessel, and Zs and A. are the atomic number and atomic weight 
of element (i). 
Die wonazation Correction 1S Given by 


_ =i 
Hp = (1- ¥(AB)) 
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In this experiment ioe was less than 1.01, i172 ere 


Landau ionization correction was less than one percent. 


D. | ERROR ANALYSIS 
ie oyStematic Errors 

Because this experiment measured the ratio of the 
Geuteron to the proton cross section, many systematic errors 
in each of these cross sections, such as solid angle, scat- 
tering angle, SEM efficiency, target thickness, and counting 
meee COrrections were effectively cancelled when the ratio of 
the two cross sections were formed. This follows because the 
data for both cross sections were taken during one continuous 
LINAC running period with the geometries the same and the beam 
Characteristics, excluding the incident energy, as nearly 
identical as possible. 

A systematic error which does not cancel in this 
ratio experiment is that which can arise because of possible 
impurities in the hydrogen and deuterium gases. The supplier 
of the target gases advertized 99.5% pure deuterium and 
99.995% pure hydrogen. However, after closer investigation 
it was. discovered that their analytic equipment can only de- 
Beer 0.3% or more contamination of D, in H, or Ho an D,. All 
the data taken in this experiment used only one bottle of 
deuterium and one of hydrogen, thus all data were taken with 
the same amount of possible contaminant. If there were a 
difference in the contamination of H, or Do this difference 


would have caused a non zero intercept of the plot of Gp 
‘ n 


2 
Versus q . For example, a 0.3% contamination of D. and a 


S/S) 


0.0% contemination of He, would have lowered the intereept 
meeabout 0.0015. 

Whenever empty target background was compared to the 
background above the peak, the deuterium background always 
lea the same spectrum distribution in both the empty and 
full cases. However, at lower a (particularly ae = 0.05) 
the hydrogen spectrum above the hydrogen peak clearly showed 
a contaminant in the gas which, if it were a low atomic 
number contaminant would have been a few percent, but if it 
were a high atomic number contaminant it would have been less 
than a few tenths of a percent. The closeness of the inter- 


cept of G versus q? to zero (we report 0.0000 + Cl COtZ tes 


En 


LF model 15) indicates that the H, and D2 contaminants were, 


2 
in fact, almost equal. The following figure shows the dif- 
ference between the Ho and D, spectra. The elastic peaks in 


miese Cases would be at about 1000 counts. 
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2.) Random Euros 

The random errors in this ratio experiment, besides 
the statistical (Poisson) counting errors, came from errors 
in determining the density of the gas in the target, vari- 
ations in the SEM efficiency during a run, energy shifts of 
Swesincident beam, leakage currents in the beam current 
integrators, and very small errors in the radiation cor- 
rections resulting from errors in determining AE. 

The statistical error for each cross section, due €o 
the EStisson distribution of the number of electrons scat- 
tered for a given number of incident electrons, was calcu- 
lated by a subroutine in the same program that calculated 
the cross section and included the statistical errors made 
in determining the background subtraction. To determine the 
@menetvonal error in the cross section caused by these sta- 
tistical fluctuations, the cross section, X, may be written 
as M 

M=Constanticy 2a. (N, =B) 
k=] k 
where 
Ny. = Total counts (count rate and voltage corrected) 
obtained at the ee Spectrometer scoceing, 

B = Average background to be subtracted (determined at 
energies above the peak) and assumed to be the same 
for each spectrometer setting, and 

M = The number of settings of the spectrometer. 


M 


Let N 2 x be the total number of counts, then: 


I 


xX = Constant ~ (N-M-B). 
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Mne variance of the cross section is 


2. (2K)* 92 4 (2X)° 92 
ON N oB B 


Q 
Il 


which gives 


Z Z eee 


O O 
- (Constant) ~ (N+M p? 


where of is the variance in the average background and oe = N. 


J 
Let Br = 1 Bs (where B. = background at ee setting) and 
a 


J is the number of spectrometer settings above the peak at 


which background was taken, 


zB. J J 
aes ae dle Ee oe 
Since B = 5 then on = & »y on. = 5 »» BS 5 i = 5. 
J 1 J J 
Thus 
of & fconshanh) < [seuc = ] 


emaeche fractional standard deviation of the cross section is 


o 1, 
ioe {ens aa ; 
7) 2 


[N-MB] 


which can also be expressed as 


o2 

oes e MB (1+ M/d) 
Vn! 2 
x oO a 


where Ny = N-MB and the first term gives the variance due to 
the true counts and the second gives the variance resulting 
from che background. z 

The errors in the measured cross sections caused by 
errors in determining the density of the gas in the target 
varied slightly from run to run, but were always less than 
about 0.1%. The pressure dropped while liquid nitrogen was 


being pumped and then rose again and stayed steady until the 
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next nitrogen fill, but the variations were almost identical 

in the hydrogen and the deuterium, and were about 0.3% of the 
measured pressure. In calcu ene average pressure at 
ger start of the nitrogen £ill or at the completion of the fill 
was the pressure actually used, depending on which was more 
Gem@steant. In calculating the gas density, of course, if the 
pressure at the start (completion) of the fill was used for 
hydrogen, then the pressure at the start (completion) of the 
fill was also used for deuterium. 

Generally, six SEM efficiencies were measured during 
a data run on each gas; one was taken at the start of the run, 
Ome at each of the four peaks, and one at the completion of 
data taking. At each point, four or more measurements were 
made, and the fractional standard deviation of a single mea- 
surement was usually less than 0.2%, giving a fractional 
standard deviation of the mean of the SEM readings of less 
than 0.08% for the six readings. 

Because of problems such as variations in line voltage, 
eres beam energy was not constant and thus the scattered energy 
also would change slightly from time to time throughout a run. 
The result of these changes was a shift of counts from one 
channel to another of the counting ladder. This was most 
noticeable when data were being taken along the steep slopes 
of elastic peaks; often a low count in one channel would be 
accompanied by a high count in an adjacent channel. The 
obvious method of reducing these shifts was to close the 
emengy defining slits more, but this increased the background 


and thus the ratio of real counts to background became 
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less - giving larger errors on the measured cross sections. 
Nevertheless, the experimenter coulda not always unambiguously 
estimate how many counts were added or subtracted from a given 
m@mestic peak due to energy shifts. 

An examination of many runs shows that there was no 
apparent energy shift problem in those runs where the cross 
section was large and the slits could be closed and still 
maintain a good peak to background ratio. In other runs 
where the slits were opened wider, it appeared that a varia- 
tion of as much as a percent could have occurred at several 
points in an elastic peak. Adjacent data points and general 
elastic peak shapes were used to smooth out these large fluc- 
tuations. It is estimated that in the poorer runs these 
energy shifts could have caused an error in the average mea- 
Sured cross section ratio of as much as 0.2%. 

There were slight leakage currents in the circuits 
used to measure and integrate the beam cureenee These cur- 
rents were reduced by improving the insulation of various 
components and by cleaning all sensitive elements. The 
leakage currents on the SEM and the faraday cup, along with 
all the circuits used to accumulate the currents, were mea- 
sured and found to be almost trivial. These measurements 
were repeated during each run and remained very constant. If 
the electron beam was perfectly constant during all parts of 
a given run, the leakage current would introduce no errors. 
Usually the electron current was held constant to within 20%; 
this would cause an error in the ratio of the cross sections 


Bt apout 0.02%. 
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Finally, there is the error in measuring the exact 
position of the peak of the elastic spectrum from the plotted 
Gata. The amount of this error varies with such things as 
beam stability, energy, and the widths of the slits but 
Seenerally was about 0.01 MeV. This error causes an error 
in the measurement of 4E which in turn causes an error in the 
Schwinger, bremsstrahlung, and ionization corrections. The 
error in the measured cross section due to this AE error for 
melethree radiation corrections was of the order of 0.05%. 

The error in AE also anace an uncertainty in the final energy 
value at which one terminates the integration of the elastic 
Peeetium.,. This error caused about a 0.03% error in the mea- 
sured cross section. 

Combining all the above errors except the statistical 
errors and the energy shift errors (Since these two vary so 
meen trom run to run), it is seen that the total fractional 
error in the measured cross section is about 0.14%. In 
Seerei1on to this, the statistical errors in the cross sections 
varied from an average of less than about 0.2% to a maximum of 
about 0.4% for a few runs, while the energy shift errors varied 
from mee aren ely zero through an average of less than 0.1% with 
about 4 runs possibly having errors as large as 0.2%. A sum- 
Mary of the errors discussed above is given in Table IV-l. 

The average error in the cross section ratio due to 
all random errors except the statistical error, Sayg Others’ 
has been calculated from the values of the average cross 
section errors in Table IV-l to be 0.24%. The total error of 


the weighted mean ratio R was then formed from 


Al 





2 Z 2 
STotai ~ “stat. * “avg.Others : (IV-1) 
TABLE IV-1 


Random Errors in Cross Section Data 


Borors Average Values of Errors 
Statistical 0.3% 

Gas Density Om L 7 

SEM Efficiency 0.08% 

Energy Shifts 0.10% 

Leakage Currents 0.014% 

AR Measurement (Radiation Corrections) O20 596 

AE Measurement (Integration) 2037. 


Two Significant sources of error were avoided in this 
work that were present in previous work [2]. The use of a gas 
target in place of a polyethylene target and the independent 
determination of the cross section ratios by each channel per- 
mitted these improvements. The gas target scatters electrons 
from the target nuclei all along the beam path, but only those 
nuclei which satisfy the eaneeeane geometry, i.e., those 
which <a scatter and recoil through both a collimator which 
is about 6 inches from the gas target and the spectrometer 
slits which are 16 inches from the target, can enter the spec- 
eeometer. Thus 1£ the experimenter first scatters from hydro- 
gen, maintains the exact geometry, and then fills the target 
With deuterium and scatters from it, he does not have to 


calculate the target thickness in order to obtain a perfectly 


a2 





accurate ratio of cross sections - if target thickness were 
the only possible error in the cross section determination. 
This is not the case if he used a solid target which is 
eainner than the collimated entrance width to the apeeeene 
meter. Stewart [2], who used a thin polyethylene target to 
measure the deuteron and proton cross sections, estimated 
the target thickness error as 0.5%. A second significant 
error which was avoided in this work was the error experi- 
enced in determining the relative efficiencies of each of 
the ten channels. This error was quite informally estab- 
lished to be 0.5% in Ref. 2. In this present work each de- 
tector scanned the complete elastic spectra and the ratio 
of the deuteron cross section to the proton cross section 
was formed separately for the individual detectors. This 
eliminated any detector efficiency errors which would have 
arisen had the method of Ref. 2 been used. 

Having established what errors were made in the 
measured cross sections, one can determine Hour eee errors 


propagated into the value of G Let the standard deviation 


E° 
n 
Of a cross section be denoted by S5 and S, for the proton and 


deuteron. Let R= gt oP where 


oo = experimental cross section of the deuteron and 
oP . experimental cross section of the proton. 
Ga is determined from 
n ‘G 
a pee es) (Iv-14) 
pees iG D . E 
n ED Cc p 
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and G_, /G is determined eiromn 
i BE 
d “p 
s 2 2 
2 
Gp oP Ko Ga is) 
d Mott mx 
rae. “Rese Oe (IV-15) 
Gr oft KP Ge (a) 
Pp Mott m Ey d 


Selving first for the fractional error in Redvers. 


ffererrors in both cross sections we have [42]: 





a aa 
em any + mete) (IV-16) 


where Sp LS the standard deviation Of Etervraric Fk. 


is determined frem 





Tae error in the ratio of 


G 
E 
BGuation IV-15, thus the standard fr&ctional deviation of the 
% 
ratio Sp / Sn , denoted by S c. /e, ' is 
p E E 
d 6p 
S 
% oe 
Sa JG ec i (IV-17) 
Ey & 
Pp 


From Equation IV-14 we determine the variance in Ga 


ia 
to be: 2 
— G, \7/G_ (aq)\ 2 
sf [| es en (N= 1U2))) 
Gp Gp Gp we 
n d p 
and using Equation IV-15, the standard deviation in the 
neutron electric form factor is 
| nee 
eee Ree Ee 
Ss. =4s,|5-——— aE ps (IV-19) 
E R\R Gg KP Da 
n Mott m 


Up to this point in the error analysis only the 


errors in one channel have been discussed, i.e., Sp is the 
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error in R due to errors in the deuteron and proton cross 
sections for any given channel. Letting R, be the ratio 
measured by channel i and Spy its standard deviation, then 
the weighted mean of all four channels, Ra! for any given 


run is obtained by calculating 











R. 
- pals 
no Se 
p = iat i (IV-20) 
im 7 1 
2 @ 
avi leaine 
dl. 


and the standard deviation of this weighted mean ratio is 


2 il 
y 32 
1=4 R. 
a 
Ri and Sp are the experimental values of the ratio 
m 


and of the standard deviation of that ratio quoted in 
Table V-I for each run and are the values used to determine 
ell derived quantities for each run in this work, e.g., 


a , S ; G2 and all others. 


En Gp da 
n 


fee REJECTION OF DATA 

Mwenty one runs are reported in this work. In each run 
the cross section ratio, R, was measured for each of the four 
channels, giving a total of 84 measurements of R. However, 
nine of these 84 measurements have not been included in the 
data analysis because of clear indicators that each had large 


systematic errors. 


The data given in Table V-I for q? = 0.05 at @ = 60° 


does not include the data of channel 5 since the "contamina- 
tion“ background subtaction method discussed in Section III-B 
mor the q" = 0.05 data did not work well for that channel, 
Energy shifts caused clear irregularities in four runs. 
For this reason the following were not included in the data 
reported in Table V-I: channel 5 in the second q? SOTO 2g 


emennel 6 in the first a? = 0.2 run, channel 7 an tre second 


2 
q 


li 


S22, 0 = 90° run, and channels 4, 5, and 6 an the second 


2 


gi 60° run. 


O22, 9 


Il 
Il 


fre inability to establish reliable backgrounds caused 
the rejection of the data of channel 4 in both q? = 0.35 


runs. 
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Ve. EXPERIMENTAL RESULTS 


A. NEUTRON CHARGE FORM FACTORS 

The data from 21 runs are reported here. Each run con- 
sisted of four measurements of R, the ratio of the elastic 
electron-deuteron cross section to the elastic electron- 
proton cross section. Of these 84 measurements of R, 9 are 
not included in the reported data for reasons given in 
Section IV-E. 

Table V-I gives the q? transferred to the deuteron, the 
scattering angle 9, the incident energies of the electron 
maegaent on the proton, E, and deuteron, Sime the ee 
mental ratio a and the standard deviation of the ratio R as 
determined by the total error (defined in Section IV-D) and 
the statistical (Poisson) error for that run. Also given are 
gquantities.which were determined by calculations. These 
Peemiiities are explained in Section II. Using the experi- 
mental results and the calculated quantities given in Table 
V-I, the ratio of the deuteron electric form factor to the 
proton electric form factor has been derived and is listed 


Bomene Light hand column. 


The neutron electric form factor G, was obtained by 


E 
n 
using Equation II-43 in the following form 
2 
1 
- m mx p al 2 
n a, (9) KO Gp (a4) EB p 


17) 





mine deuteron charge form factor Ci listed in Appendix C and 


2 
the G, 1S. and G (qa) values given in Table V-I were used 


1 
ad "Pp p- 
mon solve for Ga TINGS Va bueno Ci depends upon the deuteren 
n 
wave function as well as q?. Three sets of the 21 values of 


Ga have been calculated for Lomon-Feshbach models #1, #5, 
n 
and #15. 


The best estimate of the weror in each of these 21 values 


of G. , for each LF model, is the error obtained using Sin 


E otal 


n 
from Equation IV-l. The best estimate of the value of 


er Ga reported for each of the three models was obtained 
n 


dq 

from these sets of CE weighted by Smotal’ For comparison, 
5 G., was also calculated from the Gs data weighted by its 
dq n a 


statistical error only, and also with equal weighting given 
to each G, to provide an insight into the sensitivity of the 
n 

result to weighting of the data. 


After calculating each op and its Smotal an average SE 


amd ats S at each q? was Calculated. A plow ect rene 


Teta... 


E. (with its Stotal) 


ft the three models. For LF model 15 the value of each Ga 
n 


average G versus eG is provided for each 


and its Srotal is also given. 


For the 21 rums reported, all values of Gp were calcu- 
n 


lated from the FORTRAN programs listed in Appendix E. In 
addition, in order to check these programs in 1L5 runs, the 


values of Gr were also calculated on a Wang Model 700 cal- 
n 


culator using programs not written by the author. These 
Wang calculations used the same raw data, but used different 
calculations and independent smoothing of the elastic peaks 


to determine Ga - The values of Ga determined by each of 
n n 
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these methods showed very good agreement, well within the 
statistical errors. This agreement is convincing evidence 
that the computer programs used in the data reduction did not 
contain errors which could bias the conclusions reported in 
this paper. The ratio R reported in Table V-I is the average 
Petor those 15 runs. 


The average value of G at each ae using Lomon-Feshbach 


E 
n 


Model #1, is given in Table V-II. 


The average value of Gy at each a uSing Lomon-Feshbach 
n 
Model #5, is given in Table V-III1. 


The value of Gp calculated from L-F Model 15 and its 
n 


Sc Hor Gach of the 21 runs is given in Table V—IV. 
i 
n 


The average value of G at each a. using Lomon-Feshbach 


E 
n 


Model #15 is given in Table V-V. 


De ROSENBLUTH PLOT 

At oe = 0.2 deuteron and hydrogen data were taken at 60, 
75, 90, and 120 degrees in order that a Rosenbluth Plot could 
emeconstructed. The theoretical interpretation of the 
Rosenbluth Plot is discussed in Section II-G. Table V-VI 


gives the experimental data and the derived quantities used 


as data for the Rosenbluth plot given in Figure V-5. 


C. CHARGE RADIUS OF THE DEUTERON 

The deuteron charge rms radius was calculated by Method 
II given by Schumacher and Bethe [44]. They have provided 
the charge radius of the deuteron from several deuteron wave 


functions. L-F #1 and L-F #5 have a charge radius of 1.971 fm 


LOZ 


Average Gr 


a. 
(im ~ ) 
05 
ORO 
20 
0.25 
o- 30 


a SO 


nN 


z ; 
at Each q’ Using Lomon-Feshbach Model #1 


Ce uaa) 


029663) 
O23 967 
0.88667 
0.86247 
O28 396) 
0.81796 


TABLE V-II 


GE 
n 


O200CZ 
020056 
0.0028 
0.0070 
Oo. 0072 
020095 


S (Total) 


Sp 
n 


O2COrsS 
O2G0L 
O27 001g) 
0.0027 
0.0014 
OS 00L2 


A linear fit to the data in Table V-II gives: 


Gh = (-0.0003 + 0.0012) + (0.0245 + 0.0059) q*, (V-1) 


Nn 


while a linear fit of G, weighted by the statistical S 


E G 
n En 
only gives: 
G, = (0.0001 + 0.0010) + (0.0231 + 0.0049)q*, (v-2) 
n 
and the 21 points, if equally weighted, give 
G = (-0.0004 + 0.0016) + (0.0263 + 0.0072)q° Cc 
E 


n 


Eeelot of the data in Table V=-II is given in Figure V-i. 
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02.015 
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Average G, and s (Total) 
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(Derived from L-F Model 1) 


versus 
2 
q 
O22 O22 
2 -2 
q (fm *) 
PIqure yi. 
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L-F 


0.4 


a} 


Average G 


2 
qd 


(fm~*) 
O205 
ORs) 
O520 
O25 
O30 
O35 


Cr (LEGS) 


0.96881 
02 23906 
0.88665 
0.86244 
0.83956 
One hy 90 


TACs Voie 


Sp 
n 


CeO CZ 
050036 
0.0028 
CO) 7.0) 
0.0074 
020096 


ae bach a using Lomon-Feshbach Model #5 


Ss (Total) 


“2 
n 


O.00i8s 
O200EE 
Ue CCiEE 
020027 
0.0014 
0.0019 


A linear fit to the data in Table V-III gives; 


CG. pam. 0003 + 0.0012) + (0.0249 + 0.0059) q" (V-4) 
n 


While a linear fit of G weighted by the statistical Sa 


n En 


E 


only gives: 


Gp fee O00L + 0.0010) + (0.0234 + 0.0049) q” (V-5) 
nh 


eoemene 21 points, if equally weighted, give 


G, = (-0.0003 + 0.0016) + (0.0266 + 0.0072) q*. (v-6) 
n 


A plot of the data in Table V-III is given in Figure V-2. 
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Figure V-2 
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A plot of the data in Table V-IV is given in Figure V-3. 
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E 


TABLET Ey 


(Total) for Bach of the 21 Runs 


n 


Calculated from Lomon-Feshbach Model #15 
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0.94002 
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at Each q? Using Lomon-Feshbach Model #15 
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A linear fit to the data in Table V-V gives: 


GCG, =me-0.0000 + 0.0012) + (0.0200 + 0.0058) q* (V-7) 
n | 

Vieeesa linear fit of Gp weighted by the statistical Sg 

: Ee, 

Siey Gives: 

Gr pamvo. 0000 + 0.0010) + (0.0194 + 0.0050) q7 (V-8) 
n 

and the 21 points, if equally weighted, give: 

G, = (-0.0003 + 0.0016) + (0.0231 + 0.0072)q* (v-9) 


19 


A plot of the data in Table V-V is given in Figure V-4. 
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Figure V-5 ~ 
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while L-F #15 has r, = 1.965 fm. The experimental . Gp (0) 
dq n 
which were determined using L-F #1, 5, and 15 are plotted 


against the square of the deuteron charge radius, - 
d 
Given by Schumacher and Bethe, in Figure V-6. Schumacher 


, aS 


notes that such a plot yields an excellent straight line 
mor all models with the correct binding energy. The inter- 


section of this line with the value of a Ga (0) determined 
dg n 


by the very accurate neutron-electron scattering work yields 


rn The value of the deuteron charge rms radius was deter- 
d 


mined from Figure V-6 to be tp = 1, 96415 270.0074 eine 
d 
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VI. CONCLUSIONS 


Beom the 75 measurements of the vation cr Ene velactic 
feet ron-deuteron cross section to the elastic elec -ren—proroy 
cross section at low q? the following conclusions are drawn: 
i) the discrepancy in the neutron charge form factor 
measured in several earlier experiments [1,2,3,4] was re- 


Solved. The value of a Gp (0) was determined to be in 


da % n | 
agreement with the thermal neutron results for all three 


deuteron wave functions (L-F 1,5,15) which were used to analyze 
the data. 


mt) the value of —S, G,, (0) ranged from (0.0245 + 
dq n 


@.0059) ene when L-F 1 was used to extract Gp from the 
n 
@eueceron data, through (0.0249 + 0.0059) cme when L-F 5 was 


eed, to (0.0200 + 0.0058) ‘ana when L-F 15 was used. Thus 


dq* n 


the very accurate thermal neutron results of =, Gp 
dq n 
0.0004 eye and was determined to be the model (of the three 


L-F 15 produced a value of be Ge (0) in best agreement with 


(0)=0.0193+ 


Gomeroered) which best describes the deuteron. The method of 
Bete iiating the LF wave functions for different models and 
the more significant differences between the three models con- 
Sidered here is discussed in Appendix C. ~° 

iii) the charge radius of the deuteron was measured to 


be r, = 1.964] + 0.0074 fm. Prior to this result the most 
d 
recent measurement [21] was rp = 1.95 + 0202 in. =the svery 
d 


accurate value of the reported charge radius of the deuteron 


was possible because all data were taken at or below 


Li 


q? = 0.35 fe and thus it could all be used to determine the 


intercept and coefficient of the q? term in the equation 


2 eo we 
Sr Soe ore) ~ 6 Trms 


iv) as a check on the results of this experiment a 
Rosenbluth plot of all of the data taken at aC = 0.2 ie 
was produced in Figure V-5. The slope of this line is 
(0.0026 + 0.0012). The theoretical slope of this line is 
én om (see Section II-G). Using the "scaling law" (Equation 


M 2 
e 2 a ee 
II-31) this slope becomes 31 (4565 M ) where Go= (G, +6, )C,- 


p Pen 
At q° = 0.2 the theoretical slope is equal to 0.000822. ‘The 


experimental slope is about 1.5 9@ from this value, 


Oat 
However, Table V~-VI shows, as does Figure V-5, that the 
second 120° point has a signifiteantly different value of a 
Eiemeao the others. If this one point is dropped, a linear 
fit to the remaining 7 points gives a value of the slope of 
0.0013 + 0.0013 in good agreement with the theoretical slope. 


This shows that the one photon exchange is valid for the 


deuteron at this a 
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APPENDIX A 
DE VREES Bo Fit 


De Vries et al. [10] analysed the proton cross section 
measurements of Bumiller et al. [11] which consisted of 58 
daca points, the 114 points of Janssens et al. [10], and the 
Seeommcs Of Chen et al. [12] as well as the 71 deuteron cross 
section measurements which de Vries presented. These data 
were used to obtain the best fit to the three pole Clementel- 
Villi [40] model for the isotopic form factors given in 
Equations A-l. This de Vries fit is valid only in the low 


2 
faeeegion. 





S = ia 
Ge = 0.50 {4 Se he )} 
lt gq ie 1 2 
jee 26.6 
Vv “Ey 
G, = 0.50 t-—5 + (1+ V, ) } 
digg pau © i 
5 
Me 
S “My “M, 
Ge = 0.44 { 5 + 5 2G = se sen) 
lt gq lt g 1 2 
sets 26.6 
V 
V wy 
G) = 2.353 s+ (1 -Vy )} (A-1) 
li oe ae i 
5 
=p 


ele? 





ihe proton Lorm Lacters are 


ees V 
eee 
p 
ec V 
ee Sut Ou 
Pp 


while the neutron form factors are 


S V 
G, = Ga - G, 
1 
ie V 
om Cee 


‘ 


De Vries started with this theoretical model given by 
Equations A-]l and computed cross sections which were then 
compared with the experimental data. The free parameters 
appearing in the model were adjusted so that a minimum in 
x? was obtained. 

mre mass o£ the p meson was taken as a free parameter 
since the width of the resonance was so wide (100 Mev), but 
the seven other free parameters were reduced to six by ap- 
peeyang the constraint 


ace e ) = =2 
( 2G, 5 = @euE ~. 


dq n q =0 

Data from the several experiments mentioned above were 
combined in several different combinations. A reduced chi- 
square was calculated for each of these. “Combinations b and 
b' were recommended by the authors as the best fits. In this 
experiment the values of the free parameters given by the b' 
fit were used since Buchanan [68] has reanalyzed all known 
peOren, neutron, and Geuteron Gata and ‘concluded that 


de Vries b’' fit best describes the proton form factors. The 


ES 





b- combination includes the data of Janssens ect al mde vie. 


values of the parameters. 


oe 228 
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oe 1. bow 
Mp2 8.463 
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APPENDIX 
THE VIRIAL COEFFICIENTS OF HYDROGEN AND DEUTERIUM 


imechis experiment the target gases were) COneati come 
metal gas chamber that was in good thermal contact with a 
liquid nitrogen reservoir which was at atmospheric pressure. 


ii@ewchamber was filled from high pressure bottles of H, or 


Z 


D. through reducing and regulating valves to a gauge pressure 


2 
Sememeout 150 psig. The density of the gas in the target was 
then computed using the virial coefficients of Ho and ny at 
77. 35K (Liquid nitrogen temperature). Since the pressure was 
iomeeored continuously throughout a data run the virial ex-— 
pansion was calculated in terms of pressure coefficients. 
Using the data of Michels and Goudeket [69] and the method 


given in the Handbook of Physics [70] the gas.density was de- 


texymanea to be 


_ PM ‘ 3 
oS vs 
1+ B'P +c" P 
where 
7 Eee? 
B' = = 0.189409 x 10 > in“/lb 
a! & Wionieessco x lon. Sinn ureee 


At 150 psig hydrogen had a density of 0.0036723 Si Sie 
and deuterium 0.0073391 Tyee 

During a run the pressure was recorded and was hen used 
to calculate (in FORTRAN program TOPRADCR) the gas density 


which in turn was used in Equation IV-2 to obtain the number 


120 


of target atoms per oe and in the bremsstrahlung and Landau 
maaration correction calculations as well as touparteraliy 
determine the cross section in Equation IV-13. 

Tie American Institute of Physics Handbook, Third Edition, 
which was published coincidently with the completion of the 
Wemerng of this thesis, has virial coefficients different 
from those given in Reference 70. The data of Reference 70 
was used for all calculations in this work. If the data of 
the AIPH, Third Edition, are used the primary results of this 


work will not be changed. That is, the slope of G the 


E e 
n 
radius of the deuteron and the slope of the Rosenbluth plot 


are unchanged. There will, however, be a difference in the 
rE (O) versus a and of the Rosenbluth 


n 
ifveote This change in virial coefficients will cause the 


intercepts of both G 


intercept of G to be raised to 0.0023, which might in- 


1) 
n 
Gicate that the “contaminent"™ in the He discussed Greseceion 


IV-D is several tenths of a percent. 


a2t 


ABPENDIX C 
E@MON-FESHBACH DEUTERON ELECTRIC CHARGE STRUCTURE FACTORS 


mene deuteron, electric charge structure factor, Cre is 


defined as 
oe Vu? Ene v? (x) |g, (x ) dr 


where u(r) and v(r) are respectively the deuteron S-and D- 
state wave functions normalized so that Ch = 1 when g = O, 

Jo is the spherical Bessel function of order zero, q is the 
four-momentum transferred and r is the separation between 
the neutron and proton. In order that one might extract the 
neutron form factor from data which give the shape of the 
deuteron one needs to have Che which means that one must 
have u(r) and v(r). These are determined from the nuclear 
model used. Lomon and Feshbach derived nuclear models from 
nucleon-nucleon data by fitting a boundary condition model 
interaction determined largely by field-theoretic forms [15]. 
Lomon and Feshbach use an energy independent boundary condition 


for the deuteron wave function, 


~ i oF 
PREY, (F Sp) 


where Ly 1S the range of the potential.) Ihis condition muse 
be different in each angular momentum state [72]. E.L. Lomon 
[14] has provided values of os based upon three different 


nuclear models. These values of C_ are listed in Table C-I. 


E 


ie 


The three models which have been considered are called LF#1, 
LFH#5 and LF#15 in agreement with Ref. 15. Model L-F #1 has 
a deuteron binding energy of 2.2242 MeV and 4.57 percent 
D-state; model L-F #5 has a binding energy of 2.2245 MeV and 
5.202 percent D-state; model L-F #15 has a binding energy of 


2.2236 MeV and 7.55 percent D-state. 


AWE Gee 


LOMON-FESHBACH DEUTERON ELECTRIC 
CHARGE STRUCTURE FACTORS 


ine) C,, (LEH) C,, (LHS) C,, (LFHL5) 
omo2s 0.98414 0.98413 0.98423 
0.05 0.96881 0.96881 0.96899 
0.10 0.93967 0.93966 0.94002 
ones auoioss 0.91234 0.91285 
0.20 "0.88667 0.88665 0.88731 
Ones 0.86247 0.86244 0.86324 
0.30 0.83961 OmesosG 0.84049 
0.35 0.81796 0.81790 0.81895 
0.40 0.79741 0.79735 0.79851 
0.45 0.77789 0277780 0.77907 
0.50 0775909 0.75921 Oe e0s7 
6.55 Owy4156 0.74147 0.74291 
0.60 Gua 2462 0.72452 0.72604 
i, 3 0.54649 0.54626 0.54854 
ew 0.43328 0.43293 0.43555 
340 0032637 0.32588 0.32867 
4.0 0.25438 omesene Ove5es7 
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APPENDIX D 


THE GAS TARGET SYotEM 


A schematic diagram of the gas target system as it was 
developed for this experiment is shown in Figure D-l. A 
detailed scale drawing of the gas target assembly is given by 
Savage [71]. In this section both the best and the poorest 
techniques which were used in this experiment to qerate the 
gas target system are discussed and hopefully this will pre- 
vent others from following some unprofitable techniques. 

As the linear accelerator was being turned on for a run 
the gas fill system was connected together between valves 
V7 and V8 and valves Vl, V2 and V7 were opened. (Actually 
V7 and V8 were always open during all the runs done with the 
final technique.) The pressure gauge was checked to see that 
it read zero gauge pressure. Then valve Vl was closed and 
valves V3, VS and V@"were opened. The forepump was turned on 
to evacuate the complete gas fill system up to valves V4 and 
V6 which are parts of the gas bottles. This pulled the system 
down to about 50 microns. 

After the accelerator had been tuned LO Gene Weorrect 
energy a large (50 liter) liquid nitrogen dewar was connected 
to feed the liquid nitrogen reservoir. Temperature sensors 
were placed in this reservoir. These sensors were taped to a 
plastic rod and were arranged along this rod as shown in 


Figure D-1. Whenever the nitrogen level fell below the top 
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Figure D-I 





sensor the Cryomiser would open the valve in the compressed 
air line forcing air into the dewar and, since the dewar had 
nmprecstire seal at the top, this fomced liguid agereg nm. 116] 
the gas target reservoir. Through trial and™@errorme2t was 
found that one minute of pumping was satisfactory; if longer 
"fill" times were set on the Cryomiser it was found that the 
gas target assembly became so cold that it would harden the 
"ering seal and a serious vacuum leak would occur. A 
shorter fill time did not adequately fill the reservoir. 

The lower sensor would sound an alarm in the end station if 
the nitrogen level fell below it, and this audio alarm was 
reproduced in the control room. 

imethe earlysruns thesqas targetwwas allowed to cool for 
at least 15 minutes before the first gas was put in the target. 
In later runs, the gas was put in immediately after the liquid 
nitrogen was pumping; both methods were satisfactory. 

Pecer a vacuum of at least 50 microns EE GE on the 
thermocouple and the liquid nitrogen had started to pump, 
valve V2 was closed. A pressure of 50 microns insured that 
any contaminants in the gas fill system would add less than 
one contaminant molecule in lo? target molecules. both regu= 
lators were turned to the closed position and valves V3 and 
V5 would also be closed. Valve V4 was then opened. The 
pressure gauge marked HP, for high pressure, read between 
2000 to 200 psig depending of course on the amount of gas 
left in the bottle. The LP, low pressure, gauge should then 
read zero. If it did valve V3 was opened. Since the gas 


target had 1 mil stainless steel windows the regulator was 


Zo 





moved slowly and steadily in order that these windows would 
not be unnecessarily strained as the pressure increased from 
vacuum to about 150 psig. This pressure was read on both the 
LP gauge and on the much more accurate, larger gauge in 
Figure III~-2 simply shown with the word pressure on its face. 
The regulator maintained the system very near the set pressure 
throughout all runs although there usually was some pressure 
drop for a short time while the gas and nitrogen systems 
reached equilibrium. A remote control TV camera was panned 
from the gas target to fee large pressure gauge after the 
operator had insured that the beam was properly aligned 
through the gas target. The camera then generally remained 
on the large pressure gauge throughout the run. Since the gas 
pressure would drop whenever a nitrogen fill was progressing, 
the Operator recorded the pressure of the gas at the start 

of the fill (i.e., when the audio warning went on in the 
control room) and he would also record the pressure at the 

end of the 1 minute fill (when the audio stopped). 

After all the hydrogen data were taken valve V4 was 
closed and valves V2 and V5 were opened. Valve V2 was opened 
very aves in order that the 1 mil stainless window would not 
be unloaded too quickly as it experienced the decrease from 
150 psig. to vacuum. Valve V5 was opened-at this time to 
ensure that no hydrogen leaked through it during the hydrogen 
fill. Again a vacuum of at least 50 microns was pulled on the 
entire system up to valves V4 and V6 (the regulators did not 
provide a vacuum barrier). After the system was down to at 


least 50 microns, valves V2, V3 and V5 were closed. The 


ney 





regulator on the deuterium was checked to ensure that it 

was closed, thus there was no step-like rise in pressure on 
the gas target qm when valve V6 was opened. Again if 
the deuterium LP gauge read zero, valve V5 was opened and the 
deuterium regulator was turned slowly to reduce the strain on 
the 1 mil steel window. After reading about 150 psig. on 

the large gauge the system had some time to come to equili- 
brium as the operator re-tuned the accelerator. After checking 
the position and appearance of the beam on the target the re- 
mote TV camera was again put on the large pressure gauge and 
the gas pressure was again recorded at the start and finish 
of each nitrogen fill as the deuterium data were taken. 

After taking the deuterium data valve V6 was closed and 
valves V2 and V3 were opened and the system was again pulled 
to at least a 50 micron vacuum. Then the accelerator beam 
was sent through the empty target in order to obtain empty 
target background at several see ena settings over the 
Gemege OL data taken on that run. 

After the accelerator was secured the forepump and the 
Giyemiser were also turned off and valve Vl was opened. This 
ie dee gas system down to air and condensed some gases, e.g. 
water vapor in the gas target. The gas target was not left 
at vacuum because it was believed that the 1 mil windows 
would fail if the target chamber should have had a vacuum 
failure. Very early in these experiments valve V8 was then 
closed in order that excess water would not condense in the 
Gas target, but as the ligmid nitrogen evaporated during the 


night the 1 mil windows were blown out - presumably due to 


Ze 





hagh vapor pressure within thie Gas tergqete 1c lates 
technique was to disassemble the fill line at valve V8. This 
protected the target from a mishap as described above should 
someone have accidently closed V8 or even Vl. 

Finally the compressed air line was separated and was 
vented directly into the liquid nitrogen reservoir directly 
above the gas target. This accelerated the rate of evapora- 
tion of the liquid nitrogen and generally would have all the 
nitrogen, and more importantly, all the water which condensed 
during this evaporation removed from the reservoir within 
12 hours. It was found that if the experimenter starts to 
fill the reservoir with liquid nitrogen for another run before 
the condensed water from the last run is completely evaporated, 
then he will most probably cause the indium seal between the 
Meas target and the reservoir to fail severely -— in fact in 
two experiments the target chamber went to air and only the 
quick closing automatic valve at the ao llaneter: saved the 
electron gun. 

It should be mentioned that the gas target with valve V8 
closed will not leak any detectable amount of gas during at 
least a tO hour period. In the early part of this work valve 
Meewes closed after the initial filling of the target with 
gas, but the pressure at the end of the run generally did 
not reproduce the initial pressure - it was generally higher. 
These pressure variations were of the order of one percent - 
too much of a variation to permit the accuracy of measurement 


required for these experiments. Thus the open target method 


Zo 


described above was used with continuous monitoring of the 
Deessure throughout the run, i.e. by recording tne pressuce 


@methe start and finish of each nitrogen fil, 
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APPENDIX E 
COMPUTER PROGRAMS USED IN DATA REDUCTION 


The following FORTRAN programs were used to calculate 
quantities explained elsewhere in this work. 

TOPRADCR calculated the Schwinger, bremsstrahlung and 
Monization (Landau) radiation corrections as well as the 
density of the deuterium and hydrogen. 

TOP#FORM calculated the experimental cross section for 
each channel. The sample listed is the program for channel 
memoer 5 (i.e. 7 =5). It was also used to obtain the 
smoothed data which was plotted to define the elastic peak. 

TOPROTON calculated the proton's Mott and Rosenbluth 
eross sections of the proton, the magnetic correction term 
and the electric and magnetic form factors of the proton 
based on de Vries’ b’ fit. 

TOPDEUTR calculated the Mott cross section of the deuteron 
and the magnetic correction term which includes within it the 


GBelativistic correction to the wave functions. 
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